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ABSTRACT 


The new NRLMSISE-00 model and the associated NRLMSIS database now include the 
following data: (1) total mass density from satellite accelerometers and from orbit 
determination, including the Jacchia and Barlier data; (2) temperature from incoherent 
scatter radar, and; (3) molecular oxygen number density, [0 2 ] , from solar ultraviolet 
occultation aboard the Solar Maximum Mission (SMM). A new component, "anomalous 
oxygen," allows for appreciable 0 + and hot atomic oxygen contributions to the total mass 
density at high altitudes and applies primarily to drag estimation above 500 kpi. 

Extensive tables compare our entire database to the NRLMSISE-00, MSISE-90, and 
Jacchia-70 models for different altitude bands and levels of geomagnetic activity. We 
also investigate scientific issues related to the new data sets in the NRLMSIS database. 
Especially noteworthy is the solar activity dependence of the Jacchia data, with whic we 
investigate a large 0 + contribution to the total mass density under the combination of 
summer, low solar activity, high latitudes, and high altitudes. Under these conditions, 
except at very low solar activity, the Jacchia data and the Jacchia-70 model indeed show 
a significantly higher total mass density than does MSISE-90. However, under the 
corresponding winter conditions, the MSIS-class models represent a noticeable 
improvement relative to Jacchia-70 over a wide range of F 10 . 7 - Considering the two 
regimes together, NRLMSISE-00 achieves an improvement over both MSISE-90 and 
Jacchia-70 by incorporating advantages of each. 
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1. Introduction 

The NRLMSISE-00 empirical model of thermospheric composition and temperature is 
now available for access and use by the scientific and operational communities (Appendex, 
subsection A.4). This paper compares the new model to the standard scientific (MSISE-90 
''''[Hedin, 1991]) and operational (Jacchia-70 [Jacchia, 1970]) empirical models presently in use, 
through statistical comparison with the previous MSIS database and with the newly added data 
sets that make NRLMSISE-00 unique. We then address scientific issues related to the newly 

added data. 


This milestone fulfills our initial goal of preserving and continuing the line of MSIS-class 
models. Empirical models of the thermosphere and upper mesosphere are an indispensable tool 
used by the operational and upper atmospheric research communities for data analysis, 
initialization of detailed physics-based models, and mission and instrument design. For over a 


decade, the models of choice among scientists have been the Mass Spectrometer - Incoherent 
Scatter Radar (MSIS-class) models of upper atmospheric composition and temperature: MSIS-86 


[Hedm, 1987], which ranges upward from 90 km, and MSISE-90, which extends from the 
ground to the exobase. Some research communities continue to use the CIRA-86 climatology, 
Ihich consists of two overlapping specifications: tables generated from MSIS-86 for the 


thermosphere (altitude z > 100 km; http://nssdc.gsfc.nasa.gov/space/model/atmos/cospar2.html)) 

and tables based on averages of global data compilations for the mesosphere and below (z < 120 

km; http://nssdc.gsfc.nasa.gov/space/model/atmos/cosparl.html). Operational communities use 

the Jacchia-class models or in some cases, the extremely limited US Standard Atmosphere 

(http://nssdc.gsfc.nasa.gov/space/model/atmos/us_standard.html). The database underlying 

operational lacchia models (1970 and earlier) consists of total mass density denved from orbital 
decay of objects which flew during the 1961-70. Past MSIS-class models derived from over two- 
decades of data on composition, temperature, and total mass density, rather than from orbital 


dynamics. 

For estimating total mass density, the NRLMSISE-00 model and its underlying database 
are at least comparable to the Jacchia-class models, given the inclusion of numerous orbital drag 
(F. Barlier, private communication to A. Hedm) and accelerometer data se<.s (F. Marcos, private 
communication). A noteworthy addition to the database is the actual orbital decay data on which 
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the Jacchia models are primarily based. The model also incorporates recent data on temperature 
and molecular oxygen number density, [0 2 ], The new data sets are extensive in size, spatial 
range and the time period covered and represent significant departures from the MSIS database 
used to produce previous generations of the model. 

This upgrade is important because the MSIS and Jacchia models do not depend on 
calendar year and cannot directly track any gradual changes in the atmosphere due to 
anthropogenic or solar influences. The only way in which empirical models can maintain 
currency with the recent state of the atmosphere is by continually adding recent data to the 
database and then modifying the parameter set. In addition, instrumentation and data processing 
methods have improved and have become more diverse, potentially allowing the addition of 
higher order terms and reducing the uncertainty of the model coefficients. To accommodate new 
data, the- formulation of the model and the methodology for generating it have become more 

robust (Appendix). 

One change is particularly worthy of mention. The inclusion of drag data m a neutral 
atmospheric model has required us to account explicitly for high altitude 0 + and hot atomic 
oxygen components, which might contribute appreciably to drag under some conditions (Section 
4.1). As a result, the NRLMSIS formulation now explicitly includes a component, called 
“anomalous oxygen,” to accoufee contribution of these species to satellite drag at high 
altitudes and permits the user to compute both the “thermospheric mass density” (or “total 
neutral mass density") provided by past generations of MSIS and an "effective” mass density, 
which denotes the sum of the thermospheric mass density and the anomalous oxygen 
contribution at altitudes near the exobase. 

To introduce the new model, Section 2 describes the newly added data sets and their 
relationship to the previous database. Section 3 presents statistical comparisons of NRLMSISE- 
00, MSISE-90, and Jacchia-70 with the NRLMSIS database. Thestatistical tables significantly 
augment a previous report on drag and accelerometer data [Hedin, 1988]. Section 4 discusses 
important scientific issues related io the new data and model, and Section 5 discusses our 
conclusions regarding NRLMSISE-00 and the direction of future development. The Appendix 
summarizes the formulation, generadon, and access of the new model and addresses the use of 
orbital drag and accelerometer data to generate NRLMSISE-00. 
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2. The NRLMSIS database and model 

2.1. Relationship to past models 

NRLMSISE-00 retains the calling sequence and arguments of MSISE-90 and earlier 
MSIS models Jn the remainder of the paper, we will often exclude the designation "E-00” for 
the sake of brevity. One can usefully interpret the NRLMSIS model as a flexible, semi-empirical 
view of its extensive underlying database. That is, the model takes statistical variability into 
account while interpolating among, or extrapolating, the underlying data sets to estimate 
composition and temperature for times, geophysical conditions, and locations not covered 
specifically by the database. As with earlier versions of the MSIS-class models, the NRLMSIS 
database includes ground-, rocket-, and satellite-based measurements. The data underlying 
MSISE-90 cover the period 1965-83 and^eriveTrom incoherent scatter radar (ISR), mass 
spectrometer, solar ultraviolet (UV) occultation, pressure gauge, falling sphere, and grenade 
detonations. Until now, the database has not included either drag measurements, on which the 
Jacchia models were based, or satellite-borne accelerometer data. The new NRLMSIS upgrade 
includes those data sets. 

Both the MSIS and Jacchia models axe sensitive to the level of geomagnetic activity and 
provide an estimate of the average upper atmospheric state under geomagnetic storm conditions. 
However, at high latitudes and high geomagnetic activity, the databases are sparse, and as 
statistical averages the models do not capture the local structure and shorter time scales 
associated with any particular storm. The NRLMSISE-00 model remains a statistical average, 
but the database now' contains more data covering extremes of location and forcing. 

2.2. Expanded database 

The most important upgrade has been the addition of recent data sets and new categories 
of data to augment the NRLMSIS dat^ase^nd model: 

(a) Satellite drag [Jacchia, 1970; Barber, 1978; Hedin, 1988]: orbit determination (1961-1971); 

(b) Accelerometer (Hedin, 1988; F. Marcos, private communication); Atmospheric Explorer 
(AE-C, D, E) MESA [Champion and Marcos, 1973], Air Force SETA [Rhoden et al. f 2000], 
CACTUS [Boucfon et al„ 1979], San Marco-5 [Ardumi et al., 1997]]; 
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(c) Incoherent scatter radar (ISR) - exospheric temperature (T„): Millstone Hill (1981-97; 
[Buonsanto and Pohlman, 1998]), Arecibo (1985-95; [Melendez- Alvira et al., 1998] 

(d) ISR - Lower thermospheric temperature (T| DW ): Millstone Hill (1988-97; Goncharenko and 
Sal ah [1998], 

(e) Solar Maximum Mission (SMM) Ch density data derived from occupation of solar UV 
emissions [Aikin et al., 1993]. 

Ingesting the drag and accelerometer data on total mass density (p) promises to remove a 
postulated deficiency of MSIS for orbital tracking applications. With the inclusion ol the Jacchia 
data, the more extensive and well-documented NRLMSIS database should equal or improve the 
statistical predictions of p and of drag over those of the Jacchia models. We expect to test this 
hypothesis by applying the model to operational precision orbit determination and prediction 
[Knowles et al., 2001], An important point is that the new data on total mass density also 
influence the model coefficients for both temperature and composition. These differences should 
become apparent as NRLMSIS is compared to additional data sets. 

The incoherent scatter radar data directly influence the model temperature, which is the 
core of the MSIS formulation. Further, the new data cover an appreciable fraction of a solar 
cycle or more and are vital for both testing the existing models and producing new versions. The 
methods of orocessing ISR data have also undergone significant improvements over the last 
decade, increasing the quality of the inferred ionospheric properties [GonzaleCimd Sulzer, 1996]. 
This factor imparts high value to our new data sets. 

The Millstone Hill data on lower thermospheric temperature (T| 0 w) cover 100 km < z 
<130 km [Goncharenko and Salah, 1998]. In this regime, the neutral temperature is 
approximately equal to the ion temperature , so that extraction of the information is 
straightforward. These high quality data permit us to check and reinforce key MSIS temperature 
model parameters and wiilike important in defining the model at the mesopause. 

The SMM mission provided data on the molecular oxygen number density [Ch] over the 
altitude range 140-200 km and over a wide range of solar activity. Prior to SMM, direct 
measurements of [O 2 ] above 150 km were not available at high solar activity. The SMM 
occultation measurements suggest that dissociation may increase sufficiently ro keep this density 
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nearly constant at 200 km as solar activity increases [Kayser, 1980’, Aikin et al, 1993]. These 
data are now part of the NRLMSIS database and will be important in determining dependence on 
the solar extreme ultraviolet (EUV) flux and on magnetic activity. Aj^a result, the data should be 
particularly useful in future analysis of better EUV proxies (e.g., Lean et al. [2001]). On the 
other hand, the longstanding conflict between mass spectrometer and solar UV occuhation 
measurements of thermospheric [Cb] has had a profound effect on NRLMSIS, because the ^ 
occultation data do not follow diffusive equilibrium, contradicting mass spectrometer data [Aikin 
et al., 1993]. At the same time, the latter could be biased toward high values by recombination 
of atomic oxygen within the instrument. As an example, the two sources disagree on average by 
a factor of two or more at 200 km. The SMM data set has therefore required alterations in the 
formulation of NRLMSIS (see Appendix) and has significantly influenced the dependence of 
[Cb] estimates on Fjo.7 (Section 4.2). 

3. Statistical comparisons of models to data 

A key component of this paper is the presentation of coarse statistical measures (bias, 5, 
and standard deviation, o) to gauge the agreement of the commonly used empirical upper 
atmospheric models (MSISE-90 and, Jacchia-70) and the new model, NRLMSISE-00. with the 
NRLMSIS database. Generally, the bias is the weighted mean of the residuals [d; - m,} between 
the data set [dj] and the corresponding model estimates [m;], thereby measuring their systematic 
differences. The weighting factor for a given data point is the squared reciprocal of the attributed 
error, and angle brackets “( )” denote the weighted mean. Positive bias (P) indicates that a model 
underestimates the measured density values on average while negative bias signifies 
overestimation. The standard deviation is (((dj - irij) 2 ) - P') A and measures a model s coverage of 
the time scales and phases implicit in the data. However, (J also reflects the noise implicit in the 
data sets, so that the standard deviation can be somewhat ambiguous. When the three models are 
compared with identical data sets, the relative values of a should give an indication of the 
agreement of the various models with measured time scales and the associated phases. 

The NRLMSIS database consists of two components: the complete data sets jicquired 
from the various sources and the subset of data ■‘selected ,, to generate the model. Hedin et a I. 
[1977] described the selection process, which was designed to ensure the widest coverage of the 
hyperspace of subroutine arguments while satisfying constraints imposed by computing and 
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storage limitations. Hedin's method also avoids dominance of the mode] coefficients by only a 
few large data sets, although this can also be accomplished by proper weighting of data-model 
residuals in computing x~. In order to insure against domination of NRLMSISE-00 by the new 
data sets and to maintain consistency with past MSIS versions, we have performed a similar 
selection from the new data. 

The tables of (P, a) values include both the model generation database and the complete, 
newly added data sets. For temperature T, the cables show p = (T,(data) -Ti (model)) and C = 
(([T,<data) - Tjtmodel)] 2 ) - p 2 ) w . For species number density and total mass density, we have 
computed the bias as P = exp(log c (p; (data) /p; (model)) -1 to provide an average fractional 
density bias and have computed the standard deviation as (J = ((loge { Pi (data) /pi (model)}) - 
p 3 ) The number of tables is sufficiently large thai they have been placed in the AGU 
Electronic Dataset Archive for access by interested readers. In our discussion below, we denote 
these tables by Ai(a)-(c) to A9(a)-(c), where “A” signifies “archive” and (a)-(c) denote, 
respectively, quiet (A p < 10) and active (A ? > 50) geomagnetic conditions and the union of all 
geomagnetic conditions (“all” data). The tables abbreviate the names of the models to N00 
(NRLMSISE-00), M90 (MSISE-90), and J70 (Jacchia-70). The rabies cover total mass density 
(p), temperature, and individual species (excluding anomalous oxygen). For the tables ielatmg 
to quiet conditions and to all data, the calculations used only data points deviating from the new 
model by 15 x the associated error, or less. For active conditions, the database was sufficiently 
sparse that all high activity points were used to construct the tables. Finally, for “all” and A p > 
50, we used the 3 hr ap inputs to the NRLMSISE-00 (“N00”) and MSISE-00 (“M90”) models, 
while for quiet conditions, we used the daily A P input. 

The bias and standard deviation are coarse measures, and are useful primarily for simple 
error analysis. Statistical' differences among models often appear to be moderate and are not 
consistent across all data sources, potentially masking trends and systematic differences. Relative 
to our database, the most obvious difference; among the models involve the standard deviation 
of the model-data residuals. While a is comparable among the models regarding temperature and 
total mass density, NRLMSISE-00 shows better performance for composition, especially as 
altitude increases. The latter effect also occurs for exospheric temperature but is not pronounced. 
For Jacchia-70, the best temperature results relative to the MSIS-class models occurred for the 


8 


combination of lower altitudes and satellite-based observations. At high geomagnetic activity, 
comparisons with the data showed similar or degraded performance of all models relative to the 
low geomagnetic activity cases. 


Consideration of key drivers or variations makes the differences among models more 
apparent (e.g., see the next section on variations with the solar EUV flux, as represented by 
F10.7). In addition, comparisons of the models for selected data sets can sometimes be 
illuminating, if the user has knowledge of those particular data sets, if the data sets derive from 
the same or similar sources, or if particular altitude or geomagnetic activity ranges are 
emphasized. As an example and as an introduction to Section 4.1 on anomalous oxygen 
contributions to total mass density, Table 1 below compares the Jacchia data set to the three 
models. As indicated above, the table shovJjhat the biases and standard deviations of the three 


,'-<A 


models are comparable in magnitude; this especially true for p, across all accelerometer and drag 

A 

data sets. However, two additional, secondary features appear. First, the Jacchia model show’s a 
consistent negative bias, on average overestimating the Jacchia data. Section 4.1 shows that this 
•is likely attributable to a nonoptimal match with F10.7 variability of the data. Second, at high 
geomagnetic activity, the standard deviations of the M90 and J70 models are consistently higher 
than that of NRLMSIS. This suggests that the new model handles spatial and seasonal 
variability somewhat better than the other models at elevated geomagnetic activity. In neither 
case, however, does (3 or a allow us to make a strong inference. 


Table 1 

STATISTICAL COMPARISON OF EMPIRICAL MODELS TO JACCHIA DATA 


A? 

Altitude 

Points 

N( 

P 

)0 

<T 

M 

P 

?0 

a 

J7 

0 

a 

< 10 

200- 400 

6236 

-0.06 

0.17 

-0.06 

0.17 

-0.04 

0.17 


400- 800 

10041 

-0.07 

0.23 

-0.08 

0.26 

-0.07 

0.25 

| 

800- 1200 

5586 

0.01 

0.23 

0.03 

0.27 

-0.05 

0.23 


> 1200 

15 

0.20 

0.09 

0.27 

0.10 

-0.18 

0.05 

Ail 

200- 400 

10456 

-0.07 

0.17 

-0.06 

0.17 

-0.07 

0.19 


400- 800 

16021 

-0.08 

0.25 

-0.07 

0.27 

-0.09 

0.28 


800- 1200 

9373 

0.01 

0.24 

0.04 

0.27 

-0.07 

0.25 




> 1200 

24 

0.22 

0.12 

0.30 

0.11 

-0.20 

0.13 

>50 

200- 400 

304 

-0.05 

0.23 

-0.07 

0.23 

-0.12 

0.25 


400- 800 

441 

-0.01 

0.36 

0.01 

0.39 

-0.17 

0.42 


800-1200 

282 

0.07 

0.35 

0.05 

0.39 

-0.14 

0.39 


4. Scientific and technical issues 

4.1. .Anomalous oxygen and solar activity 

4.1.1. Background 

‘~7 N - 

' As shown by Hedin [1989], both the Jacchia model and data from the neutral mass 
spectrometer aboard Dynamics Explorer-2 (DE-2) indicated that an appreciable hot atomic 
oxygen population could be present at high latitudes and altitudes (> 600 km) during the 
..summer. The DE-2 measurements showed an elevated oxygen component under moderate to 
high solar activity, consistent with suggestions by Yee et al [1980]. that a hot oxygen geocorona 
might exist with temperature of - 4000 K. For low solar activity, at high altitudes and summer 
high latitudes, the Jacchia-70 model showed significantly higher total mass density (and helium 
concentration) than did MSIS-86. Hedin hypothesized that a hot oxygen geocorona could cause 
vthese discrepancies. 

^ «- — 

Recent analyses of ISR data by Oliver [1997] and Oliver and Schoendorf [1999] have 

argued for a small, but non-negligible, hot oxygen component (emphasizing altitudes around 400 

km), especially at night and at the sol stices during solar minimum . Schoendorf et al[2000] have 

developed model profiles of hot oxygen for use in analyzing incoherent scatter radar data, e.g. , to 

derive T„ (Section 4.3). While the ISR sites are somewhat lower in latitude than our region of 

interest (see below), a broadening of the investigations by Oliver et al. should significantly 

augment the present understanding of the hot oxygen component and could guide our future 

upgrades of NRLMSIS. 

The emphasis on hot atomic oxygen changed significantly when Keating et al. [1998] 
analyzed neutral and ion mass spectrometer measurements aboard the Midcourse Space 
Experiment (MSX). which flew in a sun-synchronous (near-polar) circular orbit at approximately 
900 km during the most recent soJar minimum. Based on comparisons of the Jacchia-70 and 
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MSIS-86 models, similar to those of Hedin^jkeating et al. showed that the measured O'* 


concentration could account for the discrepancy in the respective model estimates of [He] and 
total mass density. 




In response to these developments, the NRLMSIS model now includes an ‘'anomalous 
oxygen" (AO) component, which represents any appreciable, persistent 0* and hot 0 
populations at higher altitudes (> 500 km). The functional form of the anomalous oxygen model 
profile is similar to that of an isothermal Chapman layer, with an adjustable magnitude and scale 
height (or tem peratu re: Appendix A. 1). The data used to evaluate these parameters were simply 
th ^dcchi^ax,d Barlier(JB ) data above 600 km. At the same time, we excluded the summer JB 
data above 600 km from the data sets used to determine the He and “cold” O components of the 
model. While the winter data above 600 km are common to He, cold O, and anomalous O, our 
tests have shown small differences between NRLMSISE-00 and MSISE-90 for these data, 
indicating less influence of the additional oxygen component on the new model during winter at 
high altitudes. Surprisingly, the Jacchia -70 model appears to agree less well with the Jacchia 
data under such conditions (next subsection). 


Our anomalous oxygen data set does not include the high altitude, spin-mode DE-2 
neutral mass spectrometer data [Hedin, 1959] for several reasons: w-*- 

(1) The drag data should account for both O* ions and hot O atoms while the DE-2^ata 
account only for the neutral atoms. Therefore, the DE-2 data could bias the fit against the 0* 
component detected by the MSX analysis. 




(2) Retrieving the hot oxygen component of DE-2 data is dependent on using a model for the 
cold oxygen component. 

(3) The DE-2 data have further limitations: small number of points (425) above 600 km, high 
solar activity (average P 10 .7 above 190), and narrow temporal coverage (only fall -winter of 
1981-2). 


Fortunately, a comparison of the Jacchia-Barlier data to the DE-2 data for high latitudes 
and elevated solar activity (81-day average: (Fio,?) >150), shows good qualitative agreement, 
implying that the drag-based data set has captured the elevated neutral density implied by DE-2. 
A tantalizing result is that the fit of the NRLMSIS anomalous oxygen component to the high- 
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altitude Jacchia-Barlier total mass density data yielded an effective temperature of approximately 
4177 K ± 3 %. We do not hold this to be definitive, given the limitations of our high altitude drag 
data set, and we have therefore chosen to maintain chis parameter at 4000 K, in line with 
previously cited references on hot oxygen. Note that such a temperature would seem to be too 
high for oxygen ions to form the primary component of our anomalous oxygen model. Clearly 
we need more data on, and modeling of, neutral and ionized atomic oxygen at high altitudes to 
develop a faithful representation of the atmosphere near the exobase. 

4.1.2. Comparison of models with high altitude Jacchia data 

Comparison of models with the Jacchia data above 600 km reveal both anticipated and 
unexpected features. First, the high altitude Jacchia data support the observations of Keating et 
al. [1998], regarding a significant enhancement in total mass density over past MSIS-class 
models, for the combination of low solar activity, high altitude, and high summer latitudes. 
However, this difference decreases rapidly with increasing F10.7 and decreasing altitude, and 
surprisingly, the Jacchia-70 model significantly overestimates the observed density at very low 
F10.7. Further, as a function of Fjq .7 under the corresponding winter conditions, MSIS-class 
models generally agree better with the Jacchia data than does the Jacchia-70 model itself. As a 
result, NRJLMSISE-00 achieves improvements over both MSISE-90 and Jacchia-70, 
incorporating advantages of each. 

Figures 1 and 2 compare the (previous day) F 10.7 - dependence of the total mass density 
from four sources: (1) the Jacchia data set on total mass density (data denoted here by Pd), (2) 
the corresponding NRJLMSISE-00 model values (Pn). (3) the MSISE-90 model values (Pm), tmd 
(4) the Jacchia-70 model values (pj). The Jacchia-70. model values provide a baseline, since the 
Jacchia models are the standard of the astrodynamics community for estimating orbital drag. The 
figures depict bin-averaged differences of natural logarithms, log ; (p.N/pj) denoted by solid lines, 
log : (p M /pj) denoted by dashed lines, and log c (p d /pj), located at the centers of the ±la vertical 
bars (<y = standard deviation). The averaging bins are 10 F|o.? units, and both abscissas and 
ordinates are averages over individual bins. The number at the center of a vertical bar is the 
approximate base 2 logarithm of the number of points in the corresponding bin. The solid 
horizontal line at an ordinate of 0.0 represents the bin-averaged Jacchia-70 values. When a 
vertical bar is approximately centered on the horizontal line, the Jacchia-70 model is in good 
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agreement with the data in the given bin. Significant displacement of a vertical bar from the 
horizontal line signifies poor performance by Jacchia-70. 

Figure 1 (a) shows the solar activity dependence of p above 900 km for the combination 
of summer and high latitudes (I9[ > 45°); the MSX data of Keating et al. [1998] correspond to 
(previous-day) F 10.7 = 71, altitude z = 900 km and 9 = 80.6° (N). For low to moderate solar 
activity (75-175), the Jacchia-70 model captures the trend and magnitude of the data somewhat 
better than does NRLMSISE-00. A comparison to MSISE-90 verifies the enhancement observed 
by Keating, but this effect diminishes as F 10.7 approaches 130. Unexpectedly, as F|q ,7 decreases 
below 75, the Jacchia-70 model overestimates the total measured mass density by an increasing 
amount. In fact, below F10.7 - 75, NRLMSISE-00 appears to agree with the Jacchia data better 
than does Jacchia-70, while above - 130 the differences between the two models are relatively 
small. For the entire data subset containing 684 points, we may also compare the usual statistical 
measures, i.e., standard deviation, <j, and bias |3 = (log c (pd/p n io<id)), where “model” signifies 
NRLMSISE-00, MSISE-90, or Jacchia-70, where brackets indicate an average over the data 
subset, and where the contribution of each datum to both averages has been weighted by the 
squared reciprocal of the attributed error. Tbe values of (P, cr) are (-0.08, 0.22) for NRLMSISE- 
00 and (-0. 10, 0.24) for Jacchia-70, implying no particular advantage to either model, in spite of 
the differences in Figure 1(a). MSISE-90 has (P, cr) = (0.19, 0.29), indicating a systematically 
low average density estimate and a poorer match of the observed F10.7 dependence. 

Figure 1(b) show's the situation for Jacchia’ s data set during summer at high latitudes in 
the altitude range 600-900 km (1085 points). Comparisons with MSISE-90 show that the effect 
observed by Keating et al. has decreased considerably in both peak value and the operative range 
of F 10 . 7 . Moreover, NRLMSISE-00 appears to give better agreement with the dependence of the 
data on F]q .7 than does Jacchia-70, and agrees especially well with the data at low solar activity. 
The statistical measures axe (p, cr) = (-0.08, 0.25) for NRLMSIS versus the Jacchia-70 values of 
(-0.02, 0.31). MSISE-90 has values of (0.05, 0.34), primarily due to poorer agreement at low 
solar activity. Notice that the bias (P) values of the models are insensitive to the striking 
differences in variation with Fjo.7 and actually attribute a modest advantage to Jacchia-70, 
while a gives a coarse indication that NRLMSISE-00 does a better job matching the FI 0.7 
dependence. 
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In winter, for the same combination of latitude and altitude, the MSIS-class models 
generally outperform Jacchia-70 when compared to the Jacchia data set as a function of Fid 7 . 
Figure 2(a), for z > 900 km, shows that the Jacchia model gives somewhat better agreement at 
low F,o .7 but varies oppositely with the data as F )0 .7 increases. At F10.7 = 98, the NRLMSIS and 
MSISE-90 models do show a 45% overestimate, but this is based on only two data points. In fact 
for Fi 0 .7 in the range SO - 120, only 1 9 data points were available, making the low F10.7 range 
difficult to evaluate. Across the entire range of solar activity, the differences among models show 
up in the model biases: ({3, a) = (-0.06, 0.19) for NRLMSISE-00 and (-0.12, 0.20) for MSISE-90 
versus the Jacchia-70 values of (-0.24, 0.35). 

higure 2(b), for z = 600-900 km, shows similar but less extreme differences in bias, with 
(P, o) = (-0.14, 0. 19) for NRLMSISE-00, (-0.22, 0.21) for MSISE-90 and (-0.20, 0.20) for 
Jacchia-70. The biases are closer primarily because the number of data points decreased with 
increasing F10.7, for which the deviation of Jacchia-70 from the Jacchia data was also increasing. 

The respective {3 and a values for NRLMSISE-00 and Jacchia-70 are remarkably similar, 
especially for the cases in Figure 1, demonstrating that statistical averaging can mask qualitative 
differences. Such filtering of model estimates by averaging over one or more arguments might 
explain the comparable performance of MSISE-90 and Jacchia-70 for “special perturbations 
(SP)” orbit determination [Marcos et ah, 1998]. The SP calculation fits a detailed numerical 
propagator (including drag) to space object observations covering a fit span of several days and 
acts to filter the density model over the fit span [Neal et al., 1998], In addition, one of the fitting 
parameters is the “ballistic coefficient,” which multiplies the atmospheric density in the drag 
term; adjusting the ballistic coefficient corrects the model bias over the fit span [Marcos et al., 
1998]. On the other hand, the detailed differences among models play a large role in the 
important function of orbit prediction , for which no observations of the space object are 
available. 

4.2. Solar ultraviolet occupation vs. mass spectrometry 

As indicated above, we have included the UV occultation observations of [0 2 ] by the 
Solar Maximum Mission (SMM: [Aikin et al., 1993]) in generating the new model. We have 
also included UV occultation data derived from a second band (channel 19) on the AE-E EUVS 
instrument. Prior to NRLMSISE-00, mass spectrometer data were more numerous than were 
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solar UV occultation data and therefore had a stronger influence on the model profile of [0 2 ] in 
the lower thermosphere. Furthermore, the UV occultation data were available only at very low 
altitudes (< 150 km), where comparisons between the two data sources were less dennitive 
[Aikin’eT al., 1993]. As a consequence, the MSISE-90 profile of [0 2 ] is in approximate diffusive 
equilibrium above - 150 km [Mewr et al, 2001], whereas observations by solar UV occultation 
differ from diffusive equilibrium by an increasing amount as altitude increases within the range 
140-240 km - a longstanding controversy [Aikin et al., 1993, and references therein]. General 
circulation models and other detailed chemistiy/dynamics model; also depart from diffusive 
equilibrium in the lower and middle thermosphere [Meier et al., 2001]. 

The SMM UV data also show weaker solar activity dependence than do the mass . 
spectrometer data [Aikfn et al", 1993]. As a result of these differences, we have modified the 
parameterization of the lower thermospheric altitude profiles of 0 2 and O to allow more 
flexibility in NRLMSISE-00, as described in the Appendix. The model can now accommodate 
solar activity dependent departures from diffusive equilibrium in the lower thermosphere. Figure 
3(a) shows that the new model compromises between the two data sources in' the altitude region 
125-225 km covered by the SMM data. Above this region the NRLMSIS [0 2 ] profile 
approaches diffusive equilibrium, and below this region, the two data sources agree as the 
atmosphere approaches a fully mixed state. 

Figure 3(b) shows the dichotomy between the newly added data and the previously 
existing data, from which the MSISB-90 model (horizontal line at 0.0) was generated. The low 
altitude AE-E UV occultation data (labeled D and E) extended the solar activity dependence of 
the MSISE-90 database. These data are reasonably consistent with MSISE-90, with mass 
spectrometer data, and with low altitude, low Fjo. 7 AE-C UV occultation dara (rigure 3(a) and 
[Aikin et al., 1993]). On the other hand, the newly added SMM data and the rocket and AE-E 
UV occultation data at 150 km (labelled K, M, and F and W, respectively) are not generally 
consistent with the MSISE-90 model, although these data appear to be consistent with each 
other. Figures 3(a) and (b) show that the magnitude of the disagreement depends on both 
altitude and Fjo. 7 t indicated above. 

Figure 3(b) also verifies that the SMM data on [0 2 ] depend more weakly on F 10 j than 
does MSISE-90, which did not include those data. Figure 3(c) demonstrates that the new 
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NRLMSISE-00 model fits the [0 2 ] database far better than MSISE-90, primarily because of the 
adjustment to the profile in the region of transition from a fully mixed state to diffusive 
equilibrium. 


This leaves open the question of accuracy for both mass spectrometry and solar UV 
occupation. Before the NRLMSIS models can properly portray the altitude dependence of [OJ, 
the community must resolve this fundamental conflict of the two major data classes. Note in the 
Appendix that the [0] profile in the lower thermosphere has also changed to accommodate both 
the SMM data and the mass spectrometer data on total oxygen number density, [O] + 2[0 2 ], The 
latter data are our primary source of information to set the [O] model in NRLMSISE-OO. 
Comparison of the new model with our total oxygen data shows an improvement in standard 
deviation, but the new model is quite similar to MSISE-90 regarding total oxygen content. 

4.3 Exospheric temperature 


T^e Millstone Hill and Arecibo incoherent scatter radar data on exospheric temperature 
^ (Ta) are kgh quality and extend the NRLMSIS database well into the 1990s. These data derive 
from ntting a model of ion heat balance and chemistry to the ion temperature profile (T|(z)), 
using ISR observables and parameterized models of neutral oxygen and temperature [e.g., 


Buonsanto and Pohlman, 1998]. The retrieval of T sx from the ISR data did not include a hoi 
oxygen component [Schoendorf et aJ. 2000], The newly added Millstone Hill data, shown in 
Figure 4(a), cover the period 1981-97. The data include the June, 1991, geomagnetic storm with 
maximum {F, 0 .7 , a p } - (250, 300} (see Litvin et al. [2000] and references therein) and another 


period around October 29, 1991 with maximum {F 10 . 7 , a,,} - {270, 235}. Litvin et al. pointed 
out that during the most intense storm period of June, 1991 , molecular ions dominated the 
chemistry , requiring a modification in the algorithms used to retrieve T ;x ; a similar situation 
apparently occurred in late October. Even after this correction, however, T c>; ranged significantly 
below the predictions of MSISE-90 for both periods. In fact, given that the difference was 
greater in October, 1991, when F l0 .7 was higher (Figure 4(a), abscissa - 275), the elevated F !0 . 7 
might also be a factor in the lower value of Ta,. 

Figure 4(b) shows the variation of the new Arecibo ISR data with previous day F 10 . 7 ; 
these data also include periods of high geomagnetic ac tivity, with some daily A p values well over 
50. Interestingly, MSISE-90 provides a somewhat better fit at high solar activity (> 240) than 



does NRLMSISE-00, suggesting that the response to solar forcing might vary with latitude. 
Unfortunately relatively few data points are available at such high values of F, 0 . ? and the 
database has not supported reirievai of a meaningful ]atitude-F 10 .i coupling term. Indeed, no clear 
advantage for either model emerges from comparison with the entire Arecibo data set on upper 
thermospheric temperature, which covers the last thirty-five years. 

The findings by Millstone Hill under elevated geomagnetic activity in 1991 led naturally 
to a search for Arecibo measurements that intersected with Millstone Hill data during 3 991. Two 
such periods occurred early in the year (mid-January and mid-March). At those times F, 0 . 7 was 
high (1 80-275) while A p was low to moderate (< 35). During these periods, the temperature bias 
relative to NRLMSISE-00 had the same sign at both ISR sites — negative for mid-January with 
CFi 07 - 180-220) arid positive for mid-March when (Fjo.? - 240-275). The respective biases were 
also similar m magnitude, though 20 - 40 % less than the biases during the June and October 
storm periods. The similarity of biases at the two sites during lower geomagnetic activity points 
to a global error source in the model; an example is the use of Fjo .7 as a proxy for the solar EUV 
flu^which drives variability of the thermospheric density on time scaJes of a day or longer. Bass 
et al. [1996] and Rhoden et al, [2000], respectively, have used drag and accelerometer data to 
explore other proxies as candidates to augment or replace F J0 . 7 . On the other hand, the differing 
signs of the bias during these lower geomagnetic activity periods contrasts with the decidedly 
negative bias observed during mid- and late 1991, when maximum a p was very high. 

Most importantly, the new ISR and total mass density data, when combined with the 
previous MS IS data sets, have changed the solar activity dependence of the temperature in 
NRLMSISE-00, relative to that of MSISE-90 (and MSIS-86), especially at higher altitudes. 

Figure 5 shows the difference in mean exospheric temperature estimates produced by the models 
as a function of latitude and (Fun). The NRLMSISE-00 T ex is above that of MSISE-90 only at 
low latitudes and for moderate to low (F 10 . 7 ) and then by only a few degrees. As solar activity 
increases above moderate values, the NRLMSISE-00 value of T„ falls below that of MSISE-90 
by a steadily increasing amount, reaching - 40 K at (F J0 . 7 ) > 220 and high latitudes, 191 > 45°. 

This difference is less pronounced at lower latitudes, as we might expect from Figure 4(b). The 
mean total mass density behaves similarly to the temperature. Checks of the individual 

NRLMSISE-00 data sets on composition, temperature, and density have generally confirmed this 
behavior. 


17 


4.4. Mesosphere inversion layers 

Mesosphere inversion layers (MIL) are regions of enhanced temperature (AT - 15-50 K) 
which have been observed in the upper mesosphere and at the mesopause at low and mid- 
latitudes and primarily at night [Meriwether and Gardner, 2000], Present theory attributes the 
phenomenon to enhancement of tidaJ structure through interaction with gravity waves; 
ultimately, a more comprehensive database of 24-hour observations in the mesosphere and lower 
thermosphere will be needed to confirm and complete the theory. Quantitative analysis of the 
MIL phenomenon is well beyond the scope of the present paper and the empirical model; here 
we merely point out NRLMSISE-00 does exhibit the shape of an upper MEL near the mesopause 
ac low latitudes during the nighttime. Under these conditions, this feature is more prominent in 
the new model than in MSISE-90 and warrants mention. Figure 6(a) shows a comparison of 
MSISii-90 and NRLMSISE-00 under these conditions while 6(b) shows a similar dayside 
comparison. The structure becomes less apparent at mid-latitudes, covers a broader altitude 
range, and appear to be less prominent at midnight than do the annual mean upper MIL profiles 
shown by Meriwether and Gardner. In this region, the NRLMSIS database contains only a few 
rocket observations (~ 60 temperature values), which do not appear to have sufficient 
information to cause this behavior in the model. The figure shows primarily that the formulation 
is sufficiently flexible to capture MIL structures in data sets. Outside of the addition of new 
Millstone Hill LTCS data in the 100-130 km altitude range, the only major difference from the 
MSISE-90 model is the imposition of hydrostatic equilibrium over a wider range (80-300 km). 

The latter factor plus the lower order tides in the model apparently have acted in concert to 
produce the MIL-Iike structure. Generating realistic MIL profiles with NRLMSIS awaits 
upgrading the model with the recent, extensive database of ground- and space-based 
observations of the upper mesosphere and mesopause. 

5. NRLMSIS Model - Present and Future 

The new database underlying the NRLMSISE-00 model incorporates data on total mass 
density (orbital drag and satellite accelerometers), recent incoherent scatter radar obsen’ations 
covering more than a solar cycle, and satellite-borne FUV occupation measurements of [O 2 J 
from SMM. The model interpolates among newly added and past data sets, often incorporating 
strengths or features of each data set. As a result, the exospheric temperature in NRLMSISE-00 
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now shows somewhat weaker dependence on F10.7 relative to MSTSE-90. In the lower 
thermosphere, the model compromises between mass spectrometer and ultraviolet occultation 
data in terms of the altitude dependence of [Oi] but follows the weaker solar activity dependence 
of UV occultation data more closely. 

The incorporation of satellite-based data on total mass density has allowed the inclusion 
of a new component - “anomalous oxygen” - to correct the model estimates of total density at 
high altitudes (near the exobase). This recognizes the conclusion of Keating et ah [1998] that 0 
can dominate drag under particular conditions and, through similar analysis, the conclusion of 
Hedin [1989] that hot oxygen could be important to drag. Comparison of NRLMSIS and the 
standard operational and scientific models to the orbit-based data of Jacchia at high altitudes has 
revealed significant differences in the seasonal and solar activity dependence of the models. The 
new model appears to provide advantages over both Jacchia-70 and MSISE-90 for estimating 
total mass density. 

The broadening of the database, along with comments by users and plans to replace or 
augment the F10.7 input with a superior index of the solar chromospheric extreme ultraviolet 
(EUV) flux, have led to modifications in the model formulation: 

1) A new coupling term between F10.7 and mean F10.7 (Appendix A.l) to permit more flexibility 
in representing the dependence on solar EUV. 

2) Anomalous oxygen model in upper thermosphere to allow for increased drag under some 
conditions; 

3) Upgraded representation of [0 2 ](z) in the lower thermosphere to allow more general and 
higher altitude departures from diffusive non-equilibrium and weaker solar activity 
dependence [Aikin et al ., 1993]; 

4) New [0](z) parameterization in the lower thermosphere to compensate for changes m the 
new definition of [0 2 ](z), primarily when fitting mass spectrometer data on total oxygen 
content, [O] -r2[0 2 ]; 

5) Hydrostatic equilibrium constraint over a wider altitude range to tie the upper and lower 
atmospheric regions together self-consistently; 

6) Nonzero thermal diffusion factor for Ar. 
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Item (5) is particularly interesting, in that NRLMSISE-00 produces a temperature 
structure similar to a mesosphere inversion layer under some conditions (e.g., equatorial region 
at night, low latitudes), even though our most proximate new temperature data are above the 
mesopause (altitude range 100-130 km [Goncharenko and Salah, 1998]). We attribute this effect 
to the self-consistent imposition of hydrostatic equilibrium across the mesopause, in combination 
with the interaction among low order tides. Presumably the next generation NRLMSIS model 
will provide a more accurate climatology of the MIL when the relevant satellite- and ground- 
based data are added to the database. 

An underlying theme of this paper and of our future work is the dependence of the upper 
atmosphere on the solar EUV flux, which is the primary driver on time scales of a day or longer. 
Section 4 shows that the dependence on the F 10.7 solar EUV proxy is different for the lespective 
empirical models favored by operational and scientific communities. The new NRLMSISE-00 
model appears to incorporate advantages of both model classes and therefore helps joclose the 
gap between these models. In addition, the atmospheric calibration method of Marcos et al. ■ 
[1998] and a follow-on implementation in terms of ultraviolet remote sensing us^near-real-time 
atmospheric data to improve density estimation for the “present” epoch [Nicholas et al., 2000]. 
Ultimately, however, the operational community seeks a better predictive capability. Marcos et 
al. [1998] have shown that the most likely route to this goal is through better solar EUV inputs to 
the models, such as that by Lean* al. [2001 J. Under funding by the NASA Living With a Star 
Program, we are now pursuing this approach. 
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Appendix: NRLMSISE-00 formulation, constraints, generation, and distribution package 
A. 1. Formulation 

The Appendix of Hedin [1987] defines the thermospheric portion of the MSIS-class 
models, for which the fundamental variable is the temperature T(z). The Bates-Walker 
temperature profile variables [Walker, 1965] are the exospheric temperature, T« ; the 
temperature at 120 km, Tuol 2 nd the temperature gradient at Zib = 120 km. These variables have 

the form (e.g., for Tno) 

T )J0 = T, 20 [1 + G 1JO (L)] , 


where the overbar signifies a global and temporal mean and the function G(L) includes constant, 
spherical harmonic, and harmonic terms, some of which are coupled and whose coefficients and 
phases represent the spatial and temporal time scales inherent in the data. In addition, G(L) 
contains polynomial or exponential terms in the solar EUV proxy (F 10.7 and (Fun). the 81-day 
average) and in geomagnetic activity. For chemical species i, the Bates-Walker profile variable 

is the number density at 120 km 

n j =n i exp[G j (L)] . 


The subscript on G distinguishes among unique coefficient sets for respective thermospheric 
variables in the model. The Bates-Walker profile represents species in thermal and diffusive 
equilibrium and includes thermal diffusion. Below a species-dependent altitude in the range 160- 
450 km, the model profiles differ from diffusive equilibrium by progressively greater amounts as 
z decreases, transitioning to a fully-mixed state at a turbopause Zh - 100 km. In that region, 
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MSIS-cIass models modify the density profile due to the effects of chemistry, dynamics, and lo»s 
and flow processes. 

The NRLMSISE-00 model incorporates the following modifications of the equations in 
Hedin [1987]: 

(1) The solar EUV dependence includes a new cross term with coefficient B * 0: 

G(Solar) = A AF {1+B A(F)) + C (AF)2 + D A(F) + E (A<F»2, 
where (F 10 . 7 ) is the 81 -day, time-centered average of Fjo/j . A(F) = (Fioa) - 150, and AF = F 107 - 

(Fio.?)- 


(2) The “anomalous” oxygen model profile, [OJ(z), represents nonthermal oxygen specie* 
(e.g., O* and hoi oxygen) inherent in the Jacchia and Barlier data sets at higher altitudes (^ 500 
km) and is similar to an ionospheric Chapman layer (e.g., Cotton et al., [1994]), 


[O a ](z) - [O a ](z lb )exp<j H(Zj b fT a) [ eX ^ j H (z a ,T a ) |_^ ^ 


: — M 


J J) 


where the geopotential 
mass of atomic oxygen, 


height is q(z, z ]b ) = z- z ]b , the scale height is H(z,T) = , m is the 

, and the constants are C = 76 km, z 5 = 550 km, T a = 4000 K, z lb = 120 km, 


and [O J(z, b ) (set by data) = 6.0 x 10 4 cm*. 

(3) Equations (A20a) and (A20b) of Hedin [1987] define a lower thermospheric density 
multiplier C, in Equation (Al2a); the purpose of this factor is to simulate chemistry and dynamic 
flow effects on various species. For [O] and [OJ this factor now takes the form 

C ' _eXP {l + exp [(z-z c )/H c ]} 

where R = Rj (1 + a, A(F)) and i = O, 0 2 . The constants are R 0 = - 0.045, Ro 2 = - 0.78, 

a f = 0.029, H c (O) = - H c (O*) = 21.2 km (Note: opposite signs), and z,(0) = z.(0 2 ) = 129.4 km. 

(4) Consistent with Ban£ and Kockarts [1973], the thermal diffusion factors in NRLMSISE-00 
are a, = -0.38 (i = He, H), 0. 17 (i = Ar), 0.0 (other species). This represents a change from 
MSISE-90, w-hich had aj = -0.4 (i = He, H) and 0.0 (other species). Pavlov [1979] and (private 
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communication, 1998) suggests the following values: -0.38 (H), -0.28 (H 2 ), -0.27 (He), 0.17 
(Ar), 0.12(0?), -0.08 (0), and 0.1 (N 2 ). 

A.2. Constraints 


For altitudes 0 < z < z th ~ 120 km, the fundamental variables define nodes and gradients 
of the temperature profile, while pressure and density are defined by hydrostatic equilibrium and 
the ideal gas law [Hedin, 1991], As mentioned above, diffusive equilibrium no longer holds for 
the MSIS-class models below altitudes - 300 km. Because we fit the temperature and individual 
species separately (different coefficient sets), the MSIS-class models do not maintain hydrostatic 
equilibrium a priori. For this reason, the mode] generation process imposes an approximate 

hydrostatic equilibrium constraint in the region 80 - 300 km. This couples the lower and upper 
atmospheric regions, modifying some details of previous MSIS versions. Finally, since the new 
data ail relate to the thermosphere, NRLMSISE-00 retained the MSISE-90 coefficients below 
while constraining coefficient values In the range 72.5 - 1 10 km to give a total mass density at 
the ground in agreement with MSISE-90. 


A. 3. Model generation 

\ 

Generating a new version of the model requires calculation of optimal values for the 
1 - 2200nonzero coefficients. Even though only a subset of the MSIS database is used to evaluate 
[Hedm et al., 1977}, the number of data points is still quite sizable (- 3 x 10 5 ), rendering an 
inclusive Levenberg-Marquardt calculation [Press et a]., 1992] compute-intensive and 
cumb ersome. Because the NRLMSIS thermospheric data are separable by mass number 


(species, temperature, rota] mass density), one can partition the process into a series of separate 
Levenberg-Marquardt (LM) y~ minimization calculations for coefficient and data subsets for 


different altitude regions, magnetic activity level, and scales of variability. Each complete series 
of LM coefficient calculations (presently numbering fifty-two) constitutes one “grand” fitting 
cycle. The grand cycles repeat until the coefficient set is stable [Hedin, 1987). This approach has 
minimized memory requirements and maximized computing speed. We eliminate severe outliers 
by selecting only data points whose residuals are less than a specified multiple (6-15) of the 
observational uncertainty. 
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The new drag and accelerometer data represent a significant departure from the above 
picture. In the present case, the total mass density provided by the model is a secondary or 
inferred quantity, given by the sum of species mass densities, apparently requiring that all of the 
species coefficients vary simultaneously to fit the data, which are extensive. To avoid this and 
other procedural difficulties, one can take advantage of the fact that different thermospheric 
species dominate the mass density in different altitude regions. Specifically, for N 2 , 0 2 , and He, 
we have used the MSISH-90 model to determine the altitude ranges where the respective mass 
fractions are greater than 50 %, thereby splitting the data into subsets. We have added theie data 
subsets to the databases supporting the individual species coefficients and have combined the 
calculations of coefficients for N 2 and exospheric temperature (T„). As described in Section 4.1, 
we have extracted the high altitude (> 600 km) Jacchia and Barlier data to compute the 
coefficients associated with anomalous oxygen species and have excluded the summer-high- 
altitude Jacchia and Barlier data in determining the standard thermospheric constituents. 

A.4. Distribution Package and Access 

The present NRLMSISE-00 distribution package is an ASCII file containing the model 
source, a test driver, and the expected output of the test driver. Users may acquire the file via 
two methods: 

(1) download from our website: http://uap-www.nrl.navy.mil/models_web/msis/msis_home.htm 

(2) send e-mail to MRLMSISE-00@uaD2.nrl.naw.mil (no subject or message), which will result 
in a reply with the file as an attachment. 
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Figures 


1 . 


2 . 


Natural logarithm of total mass density (data and model values) vs. previous-day F,o.7 for the 
combination of summer, high latitudes, and high altitudes: (a) 2 > 900 km and (b) 600 km < 
< 900 km. The density values are normalized to the Jacchia-70 model and are averaged in 
bins of 10 Fjo .7 units. Vertical bars correspond to the ±lcr range of Jacchia data within each 
bin; the Jacchia-70 model values fall on the horizontal line at 0.0. NRLMSISE-00 


corresponds to the solid curve and MSISE-90 to the dashed curve. 

Same as Figure 1 for the combination of winter, high latitudes, and high altitudes: (a) z > 
900 km and (b) 600 km <z< 900 km. 


3. 


(a) Lower thermospheric profile of [Oj]: 


locr. of ratio to MSISE-90 values. Averaging 


interval: 2 km. Solid curve: NRLMSISE-00 values. Letters and symbols identify individual 

data sources (Acronyms: Atmospheric Explorer missions, AE, Solar Maximum 

SMM) - Solar Ultraviolet Occultation : (K) SMM, (A) AE-C, 100 km, (B) AE-C, 1 30 km, 

(V) AE-C, 150 km, (D) AE-E, 100 km, (E) AE-E, 130 km, (F) AE-E, 150 km, channel 19, 

(W) AE-E, 150 km, channel 06, (M) Rocket; Mas s .. Spectromste r: (C) .AE-C, (G) AE- D, (J) 

AJE-E, (+) Rocket. 

(b) Natural logarithm of lower thermospheric [0 : ] vs. mean F, 0 .7, averaged within bins of 10 
flux units. The plot shows the mean of data values normalized by MSISE-90. Vertical bars 
correspond to the ±l<r range of noimalized [OJ values within each bin. MSISE-90 


corresponds to the horizontal line at 0.0. 

(c) Same as (b), bm with data normalized to NRLMSISE-00. which corresponds to the 
horizontal line at 0.0. 

4. T: , vs. Fjo.7 (previous day), averaged over bins of 10 flux units: (a) Millstone Hill data (1981- 
97) and (b) Arecibo data (1985-95), The quantities plotted are, respectively, differences of 
data (vertical bars) and of NRLMSISE-00 model values (solid curve) from MSISE-90 
(horizontal line at 0.0). 

5. Exospheric temperature difference between NRLMSISE-00 and MSISE-90, averted over 
longitude and time, as a function of latitude and (Fjo. 7 ) ■ 
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6. Examples of temperature profiles resembling mesosphere inversion layers, produced by 
NRLMSISE-00 (solid curve) and MSISE-90 (dashed curve). Conditions: (a) local time 2100 
hr and (b) local time 1000 hr, both at 0 UT, day 90, latitude 5°. longia.de 315% Fjo. 7 = (F.o.t) 
= 150, A p = 4. 
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Introduction to tables in AGU Electronic Dataset Archive. 

Tables of Empirical Model Comparisons with the NRLMSIS Database 

The following tables compare all data sets in the NRLMSIS database to the 
corresponding values of several empirical models: NRLMSISE-00 (MOO), MSISE-90 (M90), and 
Jacchia-70 (J70). The tables provide values of the bias or mean residual ( MEAN ) and standard 
deviation (“SD”) as functions of geomagnetic activity (quiet, high, and all levels). The tables 
include both the model generation database and the complete, newly added data sets (designated 
by footnotes). For temperature T, the tables show a mean residual = (T,(data) - T ; (model)) and 
a h (([Ti(data) - ^model)}-) - <3 2 )* in units of Kelvin. For species number density and total mass 
density, we have expressed the mean residual as an average fractional density bias p = 
exp(log £ (p; (data) /p, (model)) -1 , and have computed the standard deviation as cr = 

((log e 2 {Pi (data) /p, (model)}) - (3 2 ) 

The accompanying paper, “NRLMSISE-00 Empirical Model of the Atmosphere: 
Statistical Comparisons and Scientific Issues,” (J. M. Picone, A. E. Hedm, D. P. Drob, and A. C. 
Aikin, J. Geophys. Res. xxx, 200x), and references therein provide background on the data sets. 
Accompanying the data set same is a descriptor; abbreviations include “accel” (accelerometer), 
“drag” (drag from orbit determination), “NMS” (neutral mass spectrometer), “IMS” (ion mass 
spectrometer), “ISR” (incoherent scatter radar), and "UV occ” (solar ultraviolet occupation). 
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Tabl e 1(a) 


Comparison 
total mass 


of models to drag, accelerometer, and rocket database 
density for quiet geomagnetic activity ( A P < 10 ) 


Data Set* 

ALT 

PT 5 

M 00 

MEAN 

SD 

Jacchia 

200 - 400 

623 6 

- 0.06 

0.17 

drag 

400 - 800 

10041 

1 

O 

o 

0.23 


BOO -1200 

5 S 8 S 

o 

o 

0.23 


>1200 

15 

0.20 

0.09 


M 9 0 J 7 0 


MEAN 

SD 

MEAN 

SD ; 

l 

o 

o 

o> 

0.17 

- 0.04 

0.17 

CO 

o 

o 

1 

0.26 

- 0.07 

0.25 

0.03 

0.27 

i 

o 

o 

Ul 

0.23 

0.27 

0 . 10 

- 0.18 

0 . 05 


Barlier 

120 

- 

200 

2003 

0.03 

drag 

200 

- 

400 

9343 

0.04 


400 

- 

800 

2850 

0.12 

Cactus 

200 

- 

400 

52722 

-0.05 

accel 

400 

- 

800 

52179 

-0.05 

SETA 79 

120 

- 

200 

3792 

0.C3 

accel 

200 

- 

400 

4431 

-0 . 04 

SETA 62 

120 

- 

200 

51329 

-0.02 

accel 

200 

- 

400 

54765 

0.01 

SETA 83 

120 

• 

200 

44115 

-0.03 

accel 

200 

- 

400 

53638 

-0.01 

SETA 84 

120 

- 

200 

405B5 

-0.12 

accel 

200 

- 

400 

57258 

-0.09 

San Marco -5 

120 

- 

200 

1882 

0.37 

accel 

200 

- 

400 

26070 

0.18 


400 

- 

BOO 

4192 

0.10 

A2-C MESA 

120 

_ 

200 

23145 

0.10 

accel 

200 

- 

400 

12433 

0.03 

A£-D MESA 

120 

- 

200 

15767 

0.01 

accel 

200 

- 

400 

4427 

-0.05 

AE-E MESA - 

120 

_ 

200 

25461 

0.04 

accel 

200 

- 

400 

15190 

0.07 

rocket 

30 

- 

120 

233 

-0.04 

grenade 






rocket 

50 

- 

120 

161 

D.17 

gauge 






rocket drag 

90 

- 

120 

135 

0.03 

(sphere) 






shuttle 

60 

- 

120 

61 

-0 . 05 

accel 

120 

- 

200 

30 

0.15 


0.22 

0 . 02 

0.22 

- 0.07 

0.22 

0.21 

0.04 

0.20 

0.03 

0.22 

0.32 

0-10 

0.32 

0,17 

0.32 

0.13 

- 0.06 

0.13 

- 0.09 

0.16 

0.19 

- 0.04 

0 . 19 

- 0.07 

0.22 

0.07 

0 . 06 

0.07 

- 0.03 

0 .06 

0 . 17 

- 0 . 03 

0 . 16 ' 

- 0.04 

0 . 1 ? 

0.08 

- 0.03 

0.09 

- 0.12 

o . oe 

0 . 13 

0.02 

0 . 14 

0.00 

0.14 

0.08 

- 0.10 

0.08 

- 0.17 

0.08 

0.16 

- 0.04 

0.16 

- 0.03 

0.17 

0.08 

- 0.16 

0 . 09 

- 0.24 

0 .08 

0.43 

- 0.12 

0 . 43 , 

- 0.11 

0 .44 

0.15 

0.29 

0 .16 

0.37 

0 . 15 

0.20 

0.14 

0.19 

0.19 

0.30 

0.34 

0.12 

0.33 

0.46 

0.52 

0.17 

0.12 

0.18 

0.00 

0.19 

0.30 

0.03 

0.30 

- 0.01 

0.30 

0.15 

-0 . 03 

0.16 

- 0.01 

0.15 

0.20 

- 0.05 

0.20 

- 0.04 

0.20 

0.12 

0 .01 

0.12 

- 0.03 

0.14 

CL 24 

0 . 02 

0.23 

0.03 

0.25 

0.16 

- 0.05 

0.17 

- 0.12 

C .19 

0.23 

0.16 

0.24 

0.02 

0.28 

0 . 19 

0 . 05 

. 0.17 

- 0 . 02 

0 .19 

0 .20 

0 . 00 

0.21 

- 0.11 

0.16 

0.24 

0.21 

0.22 

- 0 . 0 B 

0,23 


complete drag and accelerometer data sets used tor calculations 



Table 1 (b) . 


Comparison of models to drag, accelerometer, and rocket database on 
total mass density for high geomagnetic activity (A e 2 =0) 


Bata Set* 

Jacchia 

drag 

ALT 

200 - 400 
400 - S00 
S00 -1200 

?TS 

304 

441 

262 

MOO 

meam 

-0.05 

-0.01 

0.07 

SD 

0.23 

0.36 

0.35 

MS 0 

mean sd 

-0.07 0.23 

0.01 0.39 

0.05 0.39 

J70 

ME AM 
-0.12 
-0.17 
-0.14 

SD 
0 . 25 
C . 42 
0 .39 

Earlier 

drag 

120 - 
200 - 
400 - 

200 

400 

800 

83 

223 

105 

-0-09 

-0.11 

-0.12 

0 . 19 
0.22 
0.47 

-0. 10 
-0.14 
-0.14 

0.20 
0 . 25 
0 .49 

-0.12 

-0.13 

-0.18 

0.20 

0.23 

0.45 

Cactus 

accel 

200 - 
400 - 

400 

300 

11487 

4336 

0 . 04 
0 . 06 

D . 2 0 
0.31 

0 . 04 

o .oe 

0 . 13 
0.31 

-0.09 

-0.15 

0 . IS 
0.29 

SETA 7 9 
accel 

120 - 
200 - 

200 

400 

4 310 
43 32 

0.00 
-0 .03 

0.11 

0.14 

0 . 01 
-0.07 

0.11 
0 . 14 

0.01 

-0.03 

0 . 09 
0 , 15 

SETA 32 
accel 

12C - 
200 - 

200 

400 

21619 

39207 

-0.04 

-0.06 

0 . 12 
0.15 

-0.10 

-0.10 

0 .14 
0 .18 

-0.05 

-0.06 

C .13 
0 . 14 

SETA 53 
accel 

120 - 
200 - 

200 

400 

14664 

26103 

-0.03 

-0.01 

0 . 09 
0.13 

-0.09 
-0 . 04 

0.0B 
0 . 1.4 

-0.10 

0.00 

0 . 09 
0 . 13 

SETA 8 4 
accel 

120' - 
200 - 

200 

400 

4 34 3 
14799 

-0 . 03 
0 . 04 

0.12 
0 . 52 

-0.10 

0.01 

0.10 

0.55 

-0.17 

-0.03 

0 . 10 
0 . 52 

San Karco-5 
accel 

200 - 
400 - 

400 

BOO 

365 

16 

0.29 
-0 . 24 

0.14 

0.41 

0 . 28 
-0.18 

0.14 
0 . 41 

0.45 
-0 . 03 

0.15 
0 .41 

AE-C MESA 
accel 

120 - 
200 - 

200 

400 

4304 

1653 

0.19 
0 . 02 

0.27 

0.27 

0.21 
0 .00 

0. 2B 
0 . 96 

0.17 

-0.05 

0 . 2 S 
0 . 97 

AS-D MESA 
accel 

120 - 
200 - 

200 

400 

192 

41 

0 .25 
0.07 

0 . 13 
0.12 

0.15 
0 .01 

0 . 12 
0.11 

0.19 

0.01 

0 . 14 

0 . 13 

AE-E MESA 
accel 

120 - 
200 - 

200 

400 

473 

133 

0.10 

0.13 

0.14 
0.3 3 

0.17 
0 . 14 

0.15 

0.39 

0 . 12 
0 . 10 

0 .21 
0.36 

rocket 

grenade 

SO - 

120 

2 

-0.22 

0.01 

-0.27 

0.02 

-0.27 

0.02 

rocket 

gauge 

so - 

120 

12 

0.49 

0 . 12 

0.31 

0 .10 

0.24 

0 .10 

rocket drag 
{ sphere) 

so - 

120 

2 

0.25 

0. 01 

0.27 

0. 00 

0.20 

0 . 01 

* Conplete 

drag and ac 

celerome 

;er data 

sets 

used for 

calculati 

on 5 
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Table 1(c). Comparison 
total mass 


of models to drag, 
density across all 


accelerornster, and rocket database 
geomagnetic activity levels 


on 


Data Set T 

Jacchia 

drag 

ALT 

200 - 400 
400 - 800 
800 -1200 
>1200 

PTS 

10456 

16021 

9373 

24 

MOO 

MEAN 

-0.07 

-0.09 

0.01 

0.22 

SD 
0.17 
0.25 
0.24 
0 . 12 

M90 

MSAN SD 

-0.06 0.17 

-0.07 0.27 

0.04 0.27 

0.30 0.11 

J7 0 

MEAN 

-0.07 

-0.09 

-0.07 

-0.20 

SD 
0 . 19 
0.28 
0.25 
0 . 13 

3arlier 

drag 

120 - 200 
200 - 400 
400 - 800 

3265 
14782 
4 5 60 

0.01 

0.01 

0.07 

0.22 

0.21 

0.31 

0.00 
0.01 
0. 08 

0.23 

0.21 

0.32 

-0.08 

0.00 

0.11 

0.22 

0.22 

0.33 

Cactus 

accel 

200 - 400 
400 - 900 

52504 

55760 

-0.05 
-0 .05 

.0.14 

0.20 

-0.04 

-0.01 

0.13 

0.20 

-0.09 

-0.10 

0 .17 
0 .25 

SETA 7 9 
accel 

120 - 200 
200 - 400 

30611 

32067 

0.01 

-0.07 

0.10 

0.15 

0.04 
-0. 05 

0.09 

0.14 

0 . 00 
-0.04 

0 . 09 
0.16 

SETA 82 
accel 

120 - 200 
200 - 400 

4 B 4 4 5 
51652 

-0.04 

-0.03 

0.10 
0 . 14 

-0.06 

-0.04 

0.11 

0.15 

-0.10 

-0.04 

0.10 
0 . 14 

SETA 03 
accel 

120 - 200 
200 - 400 

57350 

54761 

-0.07 

-0.01 

0 . 0 3 
0 .16 

-0.10 

-0.04 

0.09 

0.16 

-0.15 

-0.02 

0.09 
0 . 17 

SETA 34 
accel 

120 - 2 b0 
200 - 400 

56816. 

44996 

-0.12 

-0.09 

0.09 

0.40 

-O’. 16 
-0.11 

0.09 

0.42 

-0.22 

-D.il 

0.09 

0.42 

San Marco- 3 
accel 

120 - 200 
200 - 400 
400 - 300 

2935 

41076 

7285 

0.35 
0.21 
. 0.10 

0.13 

0.13 

0.32 

0.29 

0-17 

0.13 

0.14 

0.17 

0.32 

0.34 

0.22 

0.32 

0 . 14 
0.29 
0.51 

ae-c mesa 

accel 

120 - 200 
200 - 400 

54353 

35342 

0.11 
0 . 03 

0,17 
0 .30 

0.14 
0 . 03 

0.17 

0.30 

0 . 06 
-0.01 

0.19 

0.30 

AE-D MESA 
accel 

120 - 200 
200 - 400 

2 B 4 9 9 
3170 

0 . 01 
-0.05 

0 .16 
0.20 

-0.03 

-0.07 

0.15 

0.20 

-0.01 

-0.04 

0.15 

0.20 

ae-e mesa 

accel 

120 - 200 

200 - 400 

43799 

24332 

0 . 05 
0 .09 

0.12 

0.24 

0.03 

0.05 

0.13 

0.24 

-0.02 
0 . 04 

0,15 

0.25 

rocket 

grenade 

30 - 120 

310 

-0.02 

0.16 

-0.03 

0.16 

-0.11 

0 .18 

rocket 

gauge 

BO - 120 

247 

0.21 

0.24 

0.20 

0.25 

0.07 

0.27 

rocket drag 
(sphere) 

SO - 120 

201 

0 . 05 

0.31 

0.09 

0.27 

-0.01 

0.29 

shuttle 

accel 

30 - 120 

120 - 200 

138 

52 

-0. 07 
0.07 

0.19 

0.23 

0.00 

0.12 

0.20 

0.21 

-0.11 

-0.11 

0.17 
D . 21 

* Complete 

drag and accel erome 

ter data 

sets 

used for 

c 3 1 c u 1 a 

tions- 
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Table 2 (a) . Comparison of models to temperature database for quiet levels of 
geomagnetic activity (A P < ID) 


Data Set 

ALT 


M00 

MEAN 5D 

M90 

MEAN 5D 

J7 0 

MEAN 5D 

Arecibo 

120 - 200 

20 

5.12 

78.23 

-2.87 

82.36 

9.B5 

74.93 

I SR 

200 - 400 

32 

24.09 

4 4 . 19 

27.60 

44.91 

20.63 

43 . 72 


400 - 799 

42 

72.16 

94.17 

79.40 

9 4.31 

75.17 

91.21 


T e *‘ 

110 

31.69 

41.19 

33 . 69 

39.46 

53.17 

70.50 


T«** 

186 

-1.46 

37.93 

-3 . 86 

33.90 

-16.72 

70 . 56 


T* * r r 
1 ex 

1645 

1.50 

43.67 

-1 . 84 

41.27 

-31.96 

70.78 

Hill stone 

100 - 130’ 

2 49 B 

-0.68 

45.56 

-3.34 

43.55 

4.57 

50.76 

Hill 

m ^ 
x ex 

2 87 

41,01 

50. B2 

25.41 

45.55 

55.51 

7B.il 

I SR 

rp T T 
1 BX 

232 

-3.94 

46.00 

-30.15' 

47 . 06 

-30.77 

93 . 15 


rp 1 T t 
1 cx 

3082 

3.70 

47.74 

-16.68 

50.41 

-9.50 

84.09 

St.. Santin 

Tex 

303 

6.15 

44,78 

-12 . 72 

45.41 

11.95 

30.46 

I SR . 

Ji earner ca 

T« 

33 

-36.21 

48.46 

-35.25 

48.81 

-28 . 32 

65.62 

I SR 

Malvern 

200 - 400 

94 

-16.24 

57 . 98 

-32.77 

66.34 

37 . 10 

78.07 

I SR 

AE-C NATE 

120 - 200 

229 

-39.52 

111.18 

-44.33 

113.93 

-10.61 

111 . 94 

KM 5 

200 - 400 

212 

-37.12 

105 . 39 

-41.93 

107 . 78 

-9.05 

110. 01 

AE-D NATE 

120 - 200 

155 

-28.58 

103.34 

-30.22 

106.56 

-12.58 

107 , 63 

NMS 

200 - 400 

217 

-23.24 

94 .36 

-26.24 

94.36 

7 . 69 

107. 97 

AS-S NATE 

120 - 200 

329 

-22 . 66 

50 . 93 

-29.09 

52.23 

-10.31 

57. 84 

NMS 

200 - 400 

1003 

-36.07 

69 .99 

-32.18 

70.39 

-24.35 

78.65 

Rocket 

120 - 200 

140 

-7 . 33 

52.33 

-5.79 

55.20 

11 .39 

59.69 

NMS 

200 - 400 

82 

-25.98- 

65.35 

-32.61 

69.22 

0.52 

59.23 

DE-2 WATS 

200 - 400 

1723 

-2.38 110. 07 

-21.19 109.37 

-32.02 134.02 

NMS 

400 - 800 

358 

15.43 134.32 

-4.14 133.26 

-33.12 150.76 


T Prior to 1985 (Data used to generate MSISE-90, NRLMSISE-OO) 

** 1985-95 (Data used to generate NRLMSISS-00) 

*** 1985-95 (Entire r.ew data set) 

t Prior to 1981 (Data used to generate MSISE-90, NRLMSISS-00) 

tt 1981-9*? (Data used to generate NRLMSISS- DO } 

1 1 1 1981-97 (Entire r.ew data set) 

t 1988-97 (LTCS 2-15: Concharenko and Salah [1998]) 
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Table 2(b). Comparison of models to temperature database for high geomagnetic 
activity (A P > 50) 


Data Set 

ALT 

?TS 

MEAN 

MOO 

SD 

MEAN 

M90 

SD 

MEAN 

J7 0 

SD 

Arecibo 

120 - 200 

12 

-37.51 

75.25 

-45.49 

72.28 

-86.99 

66.17 

I SR 

200 - 400 

8 

12.16 

27 . 86 

12.98 

24.93 

-22.05 

41.95 


400 - 799 

e 

1 . 55 

4 8.31 

13.96 

46.43 

-70.11 

55.65 


P * + 

1 cx 

25 

1 . 61 

46.99 

12.93 

46 . 03 

-43 . 05 

56.16 


p w * 4 
* ex 

109 

8.68 

43.36 

18.56 

43.35 

-60. 0B 

56. 61 

Millstone 

100 - 13 0 7T 

97 

-13.24 

54.54 

-19.10 

49.92 

-2.20 

62.75 

Hill 

T 1 

1 =x 

10 

33 . 34 

40.20 

44.33 

36.52 

-6.35 

57.06 

I SR 

t w T5 

- ex 

130 

3.83 

121 . 84 

-41.76 

149.71 

12.03 

140.34 


p Til 
L cx 

45 6 

0. 45 

114 . 99 

-4 6. 39 

130.91 

10.31 

125.35 

St. Sant in 
I SR 

T 

x cx 

68 

-21 . 43 

79.14 

-24.68 

80.20 

-65.08 

6B.63 

Malvern 

ISR 

200 - 400 

7 

-49. S4 

83 . 52 

-51 . 65 

84 . 04 

-21. 09 

90.30 

AS-C NATE 

120 - 200 

62 

-37. 87 

83 . 71 

-40.32 

93.71 

-7.57 

30.37 

NMS 

200 - 400 

60 

-232.70 

253 . 42 

-239. 63 

281 . 59 

-248.31 

266.39 

AE-D NATS 

120 - 200 

14 

23 . 31 

46.07 

49.94 

42.51 

49.59 

65 . 0D 

NWS 

200 - 400 

25 

-22.25 

93 . 83 

-19.14 

100.11 

118.77 

113.96 

AE-E NATE 

120 - 200 

25 

-28 . 85 

39.79 

-32 . 13 

40.37 

-37.16 

41.07 

NMS 

200 - 400 

103 

-90 . 40 

152.53 

-36.75 

155.43 

-144 .90 

157 . 54 

DE-2 WATS 

200 - 400 

337 

-457 . 37 

602 .70 

-470.71 

600.32 

-395.19 

531.19 

NN5 

400 - 300 

190 

-22 . 90 

169.04 

-37.49 

173 . 50 

43 . 21 

189.20 


* •* 

* * 

1 

t f 
t ft 


Prior to 1985 (Data used to generate KSISE-90, NRLMSIS2-00) 
1965-95 {Data used to generate NRLMSISE- 00) 

1935-95 (Entire new data set) 

Prior to 1 9 B 1 (Data used to generate M5I3S-90, NRLMSIS3-00) 
1961-37 (Data used to generate NTRLMSISE-OD) 

1981-97 (Entire new data set) 

1900-97 (LTCS 2-15: Concharenko and Salah f!398]) 
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Table 2(c). Comparison of models to temperature database across all levels o 
geomagnetic activity 


Data Set 

ALT 

PTS 

M00 

MEAN SD 

M9 0 

MEAN SD 

J7 0 

MEAN SD 

Arecibo 

120 - 200 

55 

3.28 

80.48 

-3.46 

83.96 

-9.39 

87.31 

ISR 

200 - 400 

98 

26.28 

42,01 

28.93 

41.97 

19.82 

51.86 


400 - 793 

101 

44.97 

80.91 

58.03 

79,32 

37.90 

89.82 


T«* 

172 

31.12 

46.23 

36.19 

45.19 

34.79 

79 . 02 


Tex*’ 

411 

2 . 57 

37.57 

5.12 

38.57 

-26.29 

57 ,73 


T„“* 

3892 

2. 67 

42 . 04 

4 . 40 

41.44 

-44.64 

70,79 

Millstone 

100 - .130' 

5352 

-1.54 

46.02 

-2.70 

44.57 

11.53 

52.27 

Hill 

T ex 

559 

41 .35 

5 3.15 

30.50 

54.37 

51.30 

87.22 

ISR 

T TT 

* cot 

583 

-2 . B3 

59.14 

-24.72 

67.77 

-20.90 

94 . 93 


rn Til 

L ex 

6459 

4.76 

55.74 

-11.90 

60.29 

-4.27 

83. 18 

St. Santin 

T e x 

642 

-1.48 

47 .26 

-16.83 

46. 57 

-2 . 88 

86.92 

ISR 

Jicamarca 

T ox 

52 

-26.93 

53 . 90 

-24 . S 4 

54.29 

-38.47 

80.54 

ISR 

Malvern 

200 - 400 

170 

-20.52 

71.20 

-32.21 

70.36 

22.85 

86 . 92 

ISR 

AE-C NATE 

120 - 200 

59 B 

-39.19 

101.03 

-42 . 19 

102 . £1 

-7.73 

101 .78 

NMS 

200 - 400 

571 

-72 . 39 

117.05 

-73.10 

118 . 07 

-49.73 

127 . 07 

AE-D NATE 

120 - 200 

345 

-26 . 49 

96.02 

-24 . 52 

99 .IB 

-5.34 

101.72 

NMS 

200 - 400 

437 

-31.81 

101 .96 

-23 . 17 

101.61 

10.90 

124 .55 

AE-E NATE 

120 - 200 

664 

-17 .46 

53 .78 

-22.81 

53 . 88 

-10.82 

69 .29 

NMS 

200 - 400 

2 15 B 

-33 . 18 

70. 69 

-27 . 12 

71 . 60 

-34.88 

31 . 97 

Rocket 

120 - 200 

1 55 

-6.70 

55.47 

-6.02 

60.40 

4.44 

53 .49 

NM3 

200 - 400 

111 

-17 . 05 

73 . 76 

-23 . 69 

84 . 05 

2.30 

78 . 43 

DE-2 WATS 

120 - 200 

12 

-76.32 

548. 69 

-91 .71 

549.02 

-24.64 

523.71 

NMS 

200 - 400 

4217 

-20.43 

143.93 

-31.02 

143.27 

-36.03 

165 . 94 


400 - 800 

2213 

5.37 

142 . 22 

-5 . 85 

142.15 

-19.93 

169.57 


* Prior to 1935 (Data used to generate MSISE-90, NRLMSISE-OO) 

** 1995-95 (Data used to generate NRLMSISE- 00 ) 

+ ** 1 9 S 5 - 9 5 (Entire new data set) 

; Prior to 1931 (Data used to generate MSISE-90, NRLKSISS-QO) 

7 i 1981-97 (Data used to generate NRDMSISE-OC ) 

777 1981-97 (Entire new data set) 

i 1988-97 (LTCS 2-15: Goncharenko and Salah [1995]) 
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Table 3 {a} . Comparison of models to database on N 2 concentration ^or S’ 11 -® 
geomagnetic activity (Ap S 10) 


Data Set 

AER05-A 
N MS 

ALT 

120 - 
200 - 
400 - 

200 

400 

BOO 

PTS 

109* 
1 4 01 B * 
4 53 1 * 

MOO 

KEAN 

0.C9 

0.10 

0.24 

SD 

0.16 

0.31 

0.53 

M9 0 

MEAN 

0.11 

0.06 

0.17 

SD 

0.16 

0.31 

0.53 

J7 

MEAN 

0.10 

0.22 

0.90 

0 

SD 
0 .16 
0 . 45 
0.79 

OGO-6 

NMS 

20 D 
400 

- 

400 

eoo 

69 

1849 

0.11 

0.14 

0.17 

0.25 

0.01 

0.05 

0 .16 
0.25 

-0.05 

0.49 

0.25 
0 . 53 

San Marco-3 
NMS 

120 

200 

- 

200 

400 

9 

76 

0.05 

0.14 

0.12 

0.25 

-0.01 

0.07 

0.11 

0.25 

-0.04 

0.14 

0.14 
0 . 2 B 

St. Santin 

SO 

- 

120 

272 

-0.07 

0.43 

-0.05 

0.45 

-0 . 23 

0.40 

I SR 











Arec ibo 
I SR 

SO 

- 

120 

211 

-0.17 

0.27 

-0.17 

0.27 

-0.22 

0 . 23 

AE-C NATE 

ms 

120 

200 

400 

- 

200 
4 00 
300 

10 

1726 

234 

-0.24 

0.09 

0.26 

0.12 
0.31 
0 .40 

-0.14 

0.05 

D.1B 

0.10 

0.31 

0.39 

-0.20 

0.40 

1.03 

0.16 

0.56 

0.69 

AE-C OSS 

ms 

120 

200 

- 

200 

400 

479 
94 9 

-0.06 

-0.01 

0.15 
0 .27 

-0 . 03 
-0.04 

0.15 

0.27 

0.00 

0.16 

0.19 
0 . 3 S 

ae-d OSS 
NMS 

120 

200 

- 

200 
4 00 

256 

238 

-0.03 

-0.12 

0.16 

0,36 

-0.06 

-0.17 

0. 16 
0.37 

0.01 
0 . 02 

0.18 
0 . 50 

AE-E NACZ 
NMS 

120 

200 

400 

- 

200 

400 

BOO 

510 

1306 

838 

O.DD 
0,07 
0 .09 

0.19 
0.22 
0 ,26 

-D . 01 
-0,01 
-0 . 04 

0 .16 
0.22 
0.25 

-0.03 

-0.01 

-0.14 

0.20 
0.34 
0 .36 

Rocket 
UV occ 

120 

- 

200 

1 

0 .32 

0.00 

0 . 27 

0.00 

0-13 

0.00 

Rocket 

NMS 

B0 

120 

200 

- 

120 

200 

400 

19 

225 

87 

0 . 01 
-0.03 
-0.12 

0.30 

0.30 

0.24 

-0 . 02 
-0.0B 
-0.11 

0.30 

0.30 

0.25 

-0.22 

-0.1B 

-0.12 

0.31 

0.29 

0.26 

ESRO-4 

NMS 

200 

- 

400 

1230 

-0.12 

D . 3 1 

-0.14 

0.32 

0.01 

0 . 45 

DE-2 NAGS 
NMS 

200 

400 

- 

400 

BOO 

1955 

520 

-0.33 
-0. 26 

0.23 

0.45 

-0.41 

-0.33 

0.23 
0 . 43 

-0.37 

-0.19 

0,40 

0.64 

El SCAT 

80 

_ 

120 

23 

-0.08 

0.17 

-0 . 04 

0.21 

-0.16 

0.1] 


ISR Averages 


^Entire data set 
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Table 3 (b) . 


Data Set 


AEROS-A 

NMS 

OGO-6 
N MS 

St. Santin 
I SR 

Arecibo 
I SR 

AE-C NATE 
N ms 


AE-C OSS 
NMS 

AE-D OSS 
NHS 

AE-E NACE 
NHS 


ESRO-4 

NHS 

DE-2 NACS 
NHS 

"Entire data 


Comparison of models to database on N 2 concentration for high 
geomagnetic activity (A P > 50} 


ALT 

PTS 

MOO 

M90 

J7 0 




MEAN 

SD 

MEAN 

SD 

MEAN 

SD 

200 - 

400 

3090’ 

-0.22 

0 .86 

-0.21 

0.86 

0.02 

0 . 92 

400 - 

800 

1015* 

o 

o 

1 

0.92 

-0.09 

0,.94 

0.52 

1 .30 

200 - 

400 

5 

0.09 

0.23 

0 . 04 

0.22 

-0.23 

0.33 

400 - 

900 

224 

0.51 

0 . 57 

0.39 

0 . 60 

0.79 

0.90 

80 - 

120 

34 

0.00 

0.36 

0.00 

0.35 

-0.17 

0.31 

00 

0 

1 

120 

31 

-0 .25 

0.18 

-0.26 

0.21 

i 

o 

ro 

u> 

0.19 

2 00 - 

400 

150 

0.32 

0.56 

0.28 

0 . 57 

0.69 

0.76 

400 - 

BOO 

12 

0.23 

0.78 

0.35 

0.86 

2.79 

0,95 

120 - 

200 

128 

0.07 

0.17 

0.08 

0.17 

0.17 

0 .19 

200 - 

400 

191 

0.16 

0.31 

0.12 

6.34 

0.25 

0.54 

120 - 

200 

27 

0.29 

0.10 

0.23 

0.09 

0.36 

0.16 

200 - 

400 

26 

-0 . 05 

0.23 

-0.10 

0.25 

0.52 

0.38 

120 - 

200 

26 

-0.01 

0.16 

0 .08 

0.15 

-0 . 01 

0.19 

200 - 

400 

145 

0.24 

0.38 

0.21 

0.37 

-0.13 

0.39 

400 - 

BD0 

54 

0.24 

0.20 

0.18 

0.20 

-0.28 

0.35 

200 - 

400 

231 

0.08 

0.33 ■ 

0.07 

0.36 

0.26 

0.50 

200 - 

400 

570 

-0.30 

0.36 

-0.36 

0.36 

-0.15 

0.50 

400 - 
set 

800 

253 

0 .30 

1.79 

0.10 

1.74 

1.59 

2.21 
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Table 3(c). Comparison of models to database on N 2 concentration across all 
geomagnetic activity levels 


Data Set 

ALT 

PTS 

MOO 

M90 

J70 

AZROS-A 

120 


200 

130* 

MEAN 

0.09 

SD 

0.15 

MEAN 

0.11 

SD 

0.15 

MEAN 

0,12 

SD 

0.16 

ms 

200 

- 

400 

3 83 5 2 * 

0 . 04 

0.34 

0.04 

0.34 

0.21 

0 . 49 


4 00 

- 

800 

12342* 

0.11 

0.56 

0.12 

0.56 

0.69 

0.91 

OGO-S 

200 

_ 

400 

135 

0.09 

10.19 

0.01 

0.18 

-0.11 

0.34 

NMS 

400 

- 

800 

3578 

0 . 13 

0,26 

0.08 

0.26 

0.57 

0.68 

San Marco- 3 

120 

- 

200 

16 

-0 . 02 

0.14 

-0,07 

0.14 

-0.11 

0.18 

NMS 

200 

- 

4 C 0 

121 

0.12 

0.21 

0 .06 

0.21 

0.09 

0.24 

St. Santin 

00 

_ 

120 

558 

-0.08 

0.41 

-0.05 

0.42 

-0.23 

0.38 

I SR 

Arecibo 

80 

- 

120 

403 

-0.17 

0.28 

-0.17 

0.29 

-0.23 

0.25 

ISR 

AS-C NATS 

120 

- 

200 

10 

-0.23 

0.12 

-0 .12 

0.10 

-0.20 

0.16 

NMS 

200 

- 

400 

3149 

0.10 

0.33 

0.09 

0.34 

0.43 

0 . 61 


400 

- 

BOO 

472 

0.24 

0.45 

0 . 24 

0.46 

1.20 

0.81 

AE-C OSS 

120 

- 

200 

1180 

-0.05 

0.16 

-0.02 

0.17 

0.02 

0.21 

NMS 

200 

- 

400 

2271 

0 . 01 

0.28 

0.00 

0.28 

0.17’ 

0 . 42 

AE-D OSS 

120 

- 

200 

507 

0.00 

0.16 

-0 . 04 

0.16 

0. 05 

0. IS 

NMS 

200 

- 

400 

462 

-0,13 

0.35 

-0.15 

0.37 

0. 10 

0.61 

AE-Z NACE 

120 

- 

200 

1076 

-0.01 

0,19 

-0.01 

0.17 

-0.03 

0 . 21 

NMS 

200 

- 

400 

2712 

0.07 

0.25 

0.01 

0.25 

-0.07 

0.35 


400 

- 

800 

1638 

0 . 08 

0.26 

-0.02 

0.26 

-0.17 

0.37 

Rocket 

120 

- 

200 

2 

0 . 02 

0.26 

0.02 

0.23 

-0.07 

0.19 

UV occ 

Rocket 

80 

- 

120 

32 

-0 . 05 

0.33 

-0.08 

0,35 

-0.27 

0.35 

NMS 

120 

- 

200 

297 

-0 . 05 

0.30 

-0.06 

0.29 

-0.16 

0 .29 


200 

- 

400 

116 

-0.20 

0.55 

-0.19 

0.55 

-0.21 

0 . 50 

ESRO-4 

200 

- 

400 

2621 

-0.10 

0.32 

-0.11 

0.34 

0.08 

0 .47 

NMS 

DE-2 NACS 

200 


400 

4971 

-0.37 

0.27 

-0.39 

0.26 

-0.31 

0.44 

NMS 

400 

- 

800 

1544 

-0.26 

0.47 

-0.32 

0.47 

-0.05 

0.71 

SI SCAT 

00 

_ 

120 

42 

0.00 

0.17 

0.01 

0.18 

-0.16 

0 .14 

ISR Averages 

’’Entire data 

set 
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Table 4(a). 


Data Set 

ASROS-A 

NMS 


OGO- 6 
NMS 

San Marco- 3 
NMS 

AE-C NATS 
NMS 


AS-C OSS 
NMS 

AS-D OSS 
NMS 

AE-E NACS 
NMS 


ZSRO-4 

NMS 

DE-2 NACS 
NMS 


Comparison of models to database on total oxygen concentration 
(0T2O 2 )for quiet geomagnetic activity (A p < 10} 


ALT 

PT 5 

MOO 

M 90 

J 7 0 




MEAN 

SD 

MEAN 

SD 

MEAN 

SD 

120 

- 200 

114 - 

0.24 

0.10 

0.21 

0.10 

- 0.13 

0.11 

200 

- 400 

12 02 1 * 

0.13 

0.24 

0.14 

0.25 

0.00 

0.29 

400 

- 300 

S 725 * 

0.12 

0 . 37 

0.11 

0.37 

0.16 

0.35 

200 

- 400 

67 

0.06 

0.12 

0. 05 

0.11 

0.12 

0.15 

400 

- 600 

1325 

0 . IB 

0.16 

0.17 

0.16 

0.17 

0.20 

120 

- 200 

21 

- 0.10 

0.12 

- 0.13 

0.11 

- 0.23 

0.10 

200 

- 400 

62 

- 0.04 

0.16 

- 0.06 

0 . 15 

- 0 . 09 

0.17 

200 

- 400 

1635 

-0 . 03 

0.13 

- 0.10 

0.17 

- 0.15 

0.22 

400 

- 600 

272 

- 0 . 12 

0.22 

- 0.15 

0.20 

- 0.13 

0.24 

120 

- 200 

330 

0 . 13 

0.16 

0.03 

0.16 

- 0.15 

0.24 

2 0 D 

- 400 

273 

0 . 05 

0.24 

0.04 

0.24 

- o . o a 

0.29 

120 

- 200 

206 

0.15 

0.14 

0 . 08 

0.14 

0.02 

0.22 

200 

- 400 

214 

-D . 03 

0.25 

-0 . 06 

0.26 

- 0.11 

0.31 

120 

- 2 0 D 

292 

- 0 . OS 

0.14 

- 0.12 

0.13 

- 0.26 

0.17 

200 

- 400 

1411 

- 0.10 

0.13 

- 0.11 

0.19 

- 0. 05 

0.23 

400 

- 600 

393 

- 0 . 15 

0.18 

- 0 . 17 

0.17 

- 0.03 

0.18 

200 

- 4 DO 

1263 

- 0.17 

0 . 24 

- 0.17 

C , 24 

- 0.26 

0.27 

200 

- 4 00 

2030 

- 0.12 

0.17 

- 0,12 

0.16 

- 0.09 

0.26 

400 

- 600 

13 B 1 

- 0.06 

0.2 4 

- 0.09 

0.24 

- 0.08 

0.26 


Table 4 (b) . 


Data Set 

AERO S-A 
NMS 

OGO-S 

ms 

AE-C NATE 
NMS 


AE-C OSS 
NMS 

AE'-I)' OS'S 
NMS 

A3-S NACE 
NMS 


ESRO-4 
NMS 

DE-2 NAC5 
NMS 

♦Entire dat 


, - j t v k? c p op t o c a 1 dxvcj en c oncsn l 2 t a t x on 

Comparison or models to da^anase on 

(0+20 2 ) for high geomagnetic activity levels (A p ^ 50} 


ALT 

??s 

MOO 

MEAN 

SD 

M30 

MEAN 

SD 

J70 

mean 

SD 

200 - 
400 - 

400 

BOO 

2096* 

1553* 

0 . 05 
-0.12 

0 . 49 
D .50 

0.07 

-0.08 

0.50 

0.51 

-0 .10 
-0.13 

0 . 60 

0 . 56 

200 - 
4 00 - 

400 

600 

12 

137 

-0.01 
0 . IB 

0 . 14 
0 .18 

0.01 
0 . 19 

0.16 

0.19 

0.04 

0.18 

0.25 
0 . 22 

200 - 
400 - 

400 

600 

203 

13 

-0 . 06 
-0.26 

0.32 
0 .16 

-0.05 

-D.29 

0.32 
0 . 13 

-0.19 

-0.25 

0 . 41 
0 . 17 

120 - 
200 - 

200 

400 

142 

169 

0 .23 
0.12 

0.18 

0.30 

0.15 

0.11 

0 . 15 
0.31 

-0.05 

-0.18 

0 .28 
0.48 

1 20' - 
200 - 

2 O'O' 
400 

20- 

10 

0-72"? 
0. 0B 

0 717" 

0.19 . 

0:02- 
0 . 00 

OtIB'- 

0.20 

- 0 r 05 
-0.16 

0.2-6 
0 . 3B 

120 - 
200 - 
400 - 

200 
4 0 0 
600 

7 

133 

40 

-0 . 06 
-0.04 
-0.10 

0. 15 
0 . 32 
0.20 

-0.02 

-0.04 

-0.12 

0.15 

0.32 

0.20 

o o o 

o o o 
>-* a> 

0 .27 
0 . 43 
0 . 21 

200 - 

400 

234 

-0.11 

0.23 

-0.11 

0.23 

-0 . 29 

0 . 3B 

200 - 
400 - 

400 

600 

563 

389 

-0.16 
-0 . 03 

0.23 

0.26 

-0.16 
-0. 03 

0.24 

0.27 

-0.29 

-0.25 

0 . 39 
0 .39 


a set 
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Table 4 (c) . 

Compa 

ri son 

of models to database 

on total 

oxygen 

concentration 


(Ot 

2 C> 2 )across all 

geomagnet 

ic act 

ivity levels 



Da c a Set 

ALT 

?TS 

MOD 

M9D 

J7 0 






MEW4 

SD 

MEAN 

SD 

MEAN 

SD 

AEROS-A 

120 

- 

200 

13 5* 

0.23 

0 .10 

0.21 

0.10 

-0.12 

0 . 11 

NMS 

2 D 0 

_ 

400 

3 0 9 4 1 * 

D . 08 

0.25 

D. 10 

0.25 

-0.07 

0.34 


4 0 0 

- 

800 

22474* 

0.01 

o 

OJ 

CD 

0.04 

0 . 3B 

0.03 

0 . 43 

OGO-6 

200 

_ 

400 

156 

0. 01 

0 . 14 

0.01 

0.14 

0 . 05 

0.21 

NNS 

400 

- 

600 

3391 

0.16 

0. 17 

D . 1 7 

0.17 

0.15 

0.22 

San Marco- 3 

120 


2 0C 

3 1 

-0.11 

0.16 

-0.14 

0.14 

-0.23 

0 . 13 

NMS 

200 

- 

400 

105 

-0.05 

0.14 

-0.05 

0.14 

-0.09 

0 . 16 

AE-C KATE 

200 


400 

3153 

-0.09 

0. 13 

-0. 09 

0.13 

-0.16 

0.25 

NMS 

CO 0 

- 

500 

521 

-0.11 

0.22 

-0.12 

0.21 

-0.13 

0 . 25 

AE-C OSS 

120 


200 

349 

0.14 

0,23 

0.03 

0.23 

-0.14 

0 . 29 

NMS 

200 

- 

400 

2253 

0 . 08 

0.25 

0.07 

0.25 

-0. 08 

0 .33 

AE-D OSS 

120 


200 

392 

0.15 

0 . 15 

0 . 05 

0.15 

-0.02 

0 . 23 

NMS 

200 

- 

4 DO 

394 

-0 . 03 

0.24 

-0.05 

0.24 

-0.14 

0 .32 

AE-E NACE 

120 


200 

54 5 

-0.10 

0 . 14 

-0.13 

0 . 13 

-0.25 

0 . 19 

i 

t 

2 D 0 

_ 

400 

2 851 

-0.11 

0.21 

-0.11 

0.21 

-0.04 

0 . 26 


400 

- 

500 

1730 

-0.16 

0.19 

-0.17 

0.15 

-0.07 

0 . 20 

E5RO-4 

200 

- 

4 0 0 

2592 

-0.15 

0.25 

-0 .15 

0.25 

-0.27 

0.32 

NMS 











DS-2 MACS 

200 

. 

400 

5141 

-0 ,.13 

0.17 

-0,13 

0.17 

-0 . 14 

0.30 

NMS 

400 

- 

600 

3 52 5 

-0.03 

0.23 

-0. 09 

0.23 

-0.12 

0.25 


*Entire data set 
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Table 5(a). 


Data Set 

AEROS-A 

NMS 


OGO-6 

ms 

San Marco-3 
NMS 

AE-C NATE 
NMS 

AE-C OSS 
NMS 

AE-D OSS 
NMS 

AE-E NACE 
NMS 


ESRO-4 

NMS 

D3-2 NAGS 
NMS 


♦Entire data 


Comparison of models to database 
geomagnetic activity S 10) 


ALT ATS MOO 






MEAN 

SD 

120 

- 

200 

95* 

0.71 

0.32 

200 

- 

400 

11014* 

0.24 

0 , 32 

400 

- 

8 00 

5543" 

0.13 

0.39 

200 

_ 

400 

35 

0.20. 

0.15 

400 

- 

BOO 

1994 

0.23 

0,19 

120 

_ 

200 

1 3 

0 . 45 

0.27 

200 

- 

400 

91 

0.15 

0.21 

200 

_ 

400 

1705 

-0.30 

0.17 

400 

- 

300 

234 

-0.25 

0.13 

120 


2 00 

33 3 

-0.09 

0.27 

200 

- 

400 

S29 

0.07 

0 . 22 

120 

_ 

200 

192 

-0.09 

0,27 

200 

- 

400 

227 

-0 . 23 

0.38 

120 

_ 

200 

257 

0.12 

0.22 

200 

- 

400 

1416 

-0.06 

0. IB 

400 

- 

300 

B45 

-0.08 

0.19 

200 

- 

400 

1119 

-0 . 13 

0.32 

200 

. 

400 

1371 

-0.13 

0.16 

400 

- 

300 

1344 

-0.12 

0 .16 

300 


1200 

312 

-0 . 06 

0.20 


set 


on helium concentration for quiet 
M90 J70 


MEAN 

SD 

MEAN 

SD 

0.74 

0.32 

-0.45 

0 .37 

0.22 

0.32 

-0.06 

0.90 

0.15 

0.3B 

-0.13 

0 .79 

0 .19 

0.15 

0.25 

0.3B 

0.21 

0 , 19 

0.63 

0.42 

0.44 

0.28 

0.05 

0.19 

0.14 

0.22 

-0 . 13 

0 . 33 

-0.32 

0.1B 

-0.34 

0.76 

CNJ 

rn 

o 

i 

0.18 

-0.02 

0.73 

-0.0 8' 

0.28 

-0.54 

0.79 

ST 

o 

o 

0.21 

-0.29 

0 . 92 

1 

o 

m 

KJ 

0.27 

-0,19 

1 . 11 

-0.26 

0.37 

-0.24 

.1.25 

0.10 

0.22 

i 

o 

o 

CD 

0 .50 

r* 

o 

o 

l 

0.16 

-0.21 

0.45 

-0.10 

0.20 

-0.14 

0 .41 

-0.14 

0.31 

-0.34 

1 . 06 

-0.13 

0.16 

-0.12 

0.71 

i 

o 

US. 

0.16 

-0.25 

0.56 

1 

O 

o 

0.19 

-0.28 

0.39 
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Table 5 (b) 


Comparison of models no database on helium concentration for high 
geomagnetic activity (A P > 50} 


Data Set 

ALT 

RTS 

KC0 

M9 0 


J7 0 




MEAN 

SD 

MEAN 

SD 

MEAN 

SD 

AERO S - A 

200 - 400 

2 092' 

0.31 

0.95 

0.28 

0.95 

-c.u 

1.15 

NMS 

400 - 900 

1566* 

0 .26 

1 . 00 

0.25 

1.00 

-0.18 

1 .11, 

OGO-5 

200 - 4 DO 

13 

C . 03 

0,09 

0.02 

0.12 

0.45 

0 . 41 

NMS 

400 - 900 

140 

1 . 40 

1.49 

1.41 

1.50 

1.B2 

1.39 

AE-C NATE 

200 - 400 

174 

-0.30 

0.31 

-0.32 

0.31 

-0.33 

0.30 

NMS 

400 - 800 

12 

-0.30 

0.23 

-0.36 

0.22 

0.12 

0.54 

AE-C OSS 

120 - 200 

146 

-0.06 

0.45 

-0.06 

0.46 

-0,60 

0.99 

NMS 

200 - 400 

202 

-0.06 

0.44 

-0 . 08 

0.45 

-0.35 

1 . 02 

AE-D OSS 

120 - 200 

15 

-0.08 

0.31 

-0.16 

0.30 

0 . 02 

0.59 

NMS 

200 - 4D0' 

3’3" 

- Q\ 20' 

o', ii 

-O’. 2 8 

0. 29 

-0.17 

0.65 

AE-E NACE 

120 - 200 

8 

0.17 

0.11 

0.16 

0.12 

0.79 

0.26 

NMS 

2-00 - 400 

131 

-0.03 

0.21 

-0.06 

0.20 

0.20 

0.49 


400 - BOO 

43 

-0.13 

0.12 

-0.15 

D . 13 

0.02 

0.29 

ESRO-4 

200 - 400 

211 

-D.17 

0.34 

-0.18 

0.34 

-0.40 

0.76 

NMS 









DE-2 NACS 

200 - 400 

496 

-0.16 

0.26 

-0.17 

0.26 

-0.40 

0.64 

NMS 

400 - BOO 

572 

-0.13 

0.29 

-0.12 

0.33 

-0.45 

0.61 


BOO -1200 

44 

-0 . 17 

0 . 17 

-0.19 

0.18 

-0.50 

0.41 


^Entire data set 
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Table 5(c) . 

Comparison 

geomagnetic 

of models to database 
activity levels 

on helium 

concentration 

across 

Data Set 
AEHQS-A 

nms 

ALT 

120 - 200 
200 - 400 
400 - BOO 

PTS 

113 * 
30166" 
21340- 

MOO 

MEAN 

0.72 

0.22 

0.16 

SD 

0.31 

0.33 

0.35 

M90 

MEAN 
0.7 6 
0.21 
0.14 

SD 

0.31 

0.34 

0.36 

J70 

KEAN SD 

-0.45 0.36 

-0.16 0.89 

-0.15 0.77 

OGO-6 

NMS 

200 - 
4 00 - 

400 

BOO 

68 

3532 

0.17 

0.23 

0.15 

0.20 

0.14 

0.22 

0.14 

0.19 

0 . 40 
0.57 

0.43 

0.43 

San Marco-3 

ms 

120 - 
200 - 

200 

400 

26 
14 5 

0.45 

0.13 

0.23 

0.19 

0.43 

0.11 

0.23 
0 .20 

0.04 

-0.10 

0 . 19 
0.32 

A S-C NATE 

ms 

200 - 
400 - 

400 

300 

3277 

513 

-0.30 

-0.29 

0.19 

0.19 

-0.31 

-0.32 

0.19 

0.19 

-0.35 

-0.06 

0.74 

0.71 

AZ-C OSS 
NWS 

120 - 
200 - 

200 

400 

964 

221S 

-0 . 07 
0.05 

0.32 

0.27 

-0.07 

0.03 

0.33 

0.26 

-0.57 

-0.27 

0 .81 
0 .92 

AE-D OSS 
NMS 

120 - 
200 - 

200 

400 

365 

435 

-0.08 

-0.24 

0,30 

0.35 

-0.11 

-0.27 

0.29 

0.33 

-0.17 

-0.34 

1.10 

1.29 

AE-E NACE 
NMS 

120 - 
200 - 
400 - 

200 
400 
8 00 

518 
2 El 9 
1661 

0.12 

-0.05 

-0.03 

0.25 

0.17 

0.17 

0.10 
-0 . 0 6 
-0.09 

0.25 
0.17 
0.17 ' 

-0.03 

-0.14 

-0.09 

0.52 

0.47 

0.39 

ESRO-4 

NMS 

200 - 

40C 

2449 

-0 . 14 

0.31 

-0.16 

0.31 

-0.42 

1.03 

DE-2 NACS 
NMS 

200 - 
400 - 
600 - 

400 

800 

1200 

4351 

4620 

673 

-0.12 

-0.13 

-0.09 

0.18 

0.19 

0.10 

-0.12 
-0.13 
-0 /II 

0 . ia 
0.20 
0.13 

-0.18 

-0.32 

-0.35 

0.67 

0.56 

0.40 


^Entire data set 
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Table 6(a) . Comparison of models to database on argon concentration for quiet 
geomagnetic activity (A P < 10) 


Data Set 


ALT 

PTS 

MOO 

M9 0 

J70 





MEAN 

SD 

MEAN 

SD 

MEAN 

SD 

ASROS - A 

120 

- 200 

109* 

0 . 47 

0.22 

0.45 

0.23 

0.23 

0.19 

ms 

200 

- 400 

8566* 

0.36 

0.48 

0.32 

0.48 

0.27 

0.73 

San Marco-3 

120 

- 200 

11 

0.55 

0.21 

0 .50 

0.20 

-0.10 

0.23 

ms 

200 

- 4 OD 

6 6 

1.35 

0.29 

1.15 

0.28 

0.54 

0.36 

AE-C NATE 

200 

- 400 

341 

0.36 

0.40 

0.33 

0.42 

0.52 

0.79 

NMS 










AE-C OSS 

120 

- 200 

299 

-0.10 

0 .35 

-0.07 

0.38 

-0.20 

0.47 

L/J 

s 

200 

- 400 

3 

6.90 

0.56 

5.53 

0.59 

13.55 

0.59 

AE-D OSS 

120 

- 200 

159 

0.69 

0 . 81 

0.72 

0.82 

0.45 

0.79 

NMS 










AE-S NACE 

120 

- 200 

313 

-0.03 

o.ie 

*— i 

o 

o 

> 

0.18 

-0.25 

0.24 

NMS 

200 

- 400 

336 

0.15 

0 .26 

0.10 

0.25 

-0.13 

0.31 

Rocket 

120 

- 200 

111 

-0.13 

0.56 

-0.15 

0 . 64 

-0,42 

0 . 57 

NMS 










ESRO-4 

200 

- 400 

1033 

-0.15 

0 .40 

-o.ie 

0.40 

-0.12 

0.78 

NMS 










DE-2 NACS 

200 

- 4 OC 

1029 

-0.01 

0.43 

-0.16 

0 . 44 

-0.11 

0 .78 

NMS 

400 

- 800 

15 

5.00 

1 . OS 

3.70 

1-10 

9.19 

0.99 


* Entire data set 
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Table 6(b). Comparison of models tc database on argon concentr ation for high 
geomagnetic activity levels (A P > 50} 

Data Set ALT ?T$ MOO M90 J70 


MEAN 


AERO S- A 
N MS 

200 - 

4 00 

1723* 

0.58 

AS-C NATS 
NMS 

200 - 

400 

47 

0 . 54 

AS-C OSS 

120 - 

200 

112 

-0.01 

NMS 

200 - 

400 

4 

74 . $5 

AS-D OSS 
NK5 

120 - 

200 

11 

4 .51 

AE-E NACE 

120 - 

200 

13 

-0 .16 

NWS 

200 - 

400 

59 

0.55 

E5RO-4 

NMS 

200 - ’ 

400 

223 

-0.06 

DS-2 NACS 

200 - 

400 

465 

0.12 

NMS 

400 - 

BOO 

SI 

1.13 

* Entire data 

set 





SD 

1.53 

MEAN 

0.57 

SD 

1.52 

MEAN 

1.12 

SD 
1 . 81 

0.41 

0.56 

0.42 

0.12 

0 . 67 

0.50 

2.16 

0.14 

69.97 

0.50 
2 . 01 

0.43 

B8.57 

0 . 64 
1.91 

1.25 

4.27 

1 .25 

3.27 

1.25 

0.13 
0 . 79 

-0.15 

0.5S 

0 . 12 
0.76 

-0.46 

-0.07 

0.21 
0 . 69 

0.60 

-0.07 

0.59 

0.55 

1 . 07 

0.73 
1 . 00 

-0.09 

0.45 

0.76 

1.10 

1.27 

12.77 

1 . 12 
1.22 
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Table 5(c). Comparison of models to database on argon concentration across all 
geomagnetic activity levels 


Data Set 

AIT 

?TS 

MOO 


M9 0 

J7 0 




MEAN 

SD 

MEAN 

SD 

MEAN 

SD 

ASROS-A 

120 - 200 

13 0* 

0.43 

0,23 

0.45 

0.23 

0.24 

0.20 

ms 

200 - 40.0 

2253 3* 

0.29 

0.50 

0.28 

0,50 

0.40 

0 . 87 

San Marco- 

3 12 0 - 2 DO 

22 

0.53 

0.19 

0.51 

0.18 

-0.04 

0.23 

ms 

200 - 400 

109 

1.30 

0.27 

1.15 

0.27 

0.49 

0.34 

AE-C NATE 

200 - 400 

89 6 

0 . 51 

0.37 

0.51 

0.39 

0.78 

0 . 95 

NM5 









AE-C OSS 

120 - 200 

802 

-0.03 

0.43 

0. 00 

0.44 

-0.02 

0 . 39 

ms 

200 - 400 

25 

2 . 63 

2 .10 

2.57 

2 . 14 

3.37 

2 .46 

AS -D OSS 

120 - 200 

287 

0.56 

0.77 

0.61 

0 . 78 

0.40 

0,79 

ms 









AE-E NACE 

120 - 200 

54 6 

-0 . 0B 

0.25 

-0.05 

0.25 

-0.20 

0 . 29 

NKS 

200 - 400 

719 

0.16 

0 .39 

0.14 

0.39 

-0.13 

0.46 

Rocket 

SO - 120 

23 

-0.35 

0.54 . 

-0.36 

0.53 

-0.50 

0 . 52, 

NM5 

120 - 200 

168 

-0.14 

0.52 

-0.14 

0.49 

-0.36 

0 . 55 


200 - 400 

10 

-0.34 

0 . 02 

-0.27. 

0.03 

0. 05 

0 . 02 

ESRO-4 

20 D - 4 0 0 

2351 

-0.19 

0 . 44 

-0.19 

0.44 

0.00 

0 . 39 

NMS 









DE-2 NAGS 

200 - 400 

3045 

-0.05 

0 . 52 

-0.15 

0.53 

0.23 

0 .96 

NMS 

400 - aoo 

170 

1 . 22 

1.01 

0.64 

1.10 

8.32 

1.17 


’‘Entire data set 
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Table 7(a). Comparison of models to database on 0 2 concentration for quiet 
geomagnetic activity ( A r < 10) 

Data Set ALT PTS MOO K90 J70 


MEAN 


SMK 1 

120 - 200 

42794* 

-0.19 

UV Occ 

120 - 200 

913 

-0.03 


2 0 C - 4 00 

13043* 

-0.23 

AS-C OSS 

120 - 2 DO 

50 

-0.02 

ms 

200 - 400 

23 

0.07 

AE-D OSS 

120 - 200 

219 

D.24 

ms 

200 - 400 

71 

0. 11 

AE-E OSS 

120 - 200 

112 

0.2 4 

NMS 

200 - 400 

11 

0.21 

Rocket 

120 - 200 

11 

-0.12 

Absorption 




Rocket 

120 - 200 

211 

0.12 

ms 

200 - 400 

29 

-0.13 

AE-C EUVS 

100 

47 

-0.09 

UV OCC 

130 

20 

0.26 


150 

5 1 

0.19 

AE-E EUVS 

100 

114 

0 . 03 

UV Occ 

130 

77 

0 , 01 


150* 

72 

-0.21 


15 D 3 

152 

-0 . 06 


SD 

MEAN 

SD 

HE .AN 

SD 

0.37 

-0.40 

0.37 

-0.63 

0.41 

0.25 

-0.37 

0.29 

-0.66 

0.35 

0,62 

-0.47 

0.71 

-0.69 

0.78 

0 . 22 

-0.10 

0.24 

-0.29 

0.23 

0 . 42 

0.00 

0.39 

-0.12 

0.43 

0.34 

0.19 

0.34 

0.01 

0.36 

C .36 

0.13 

0.36 

0.05 

0.42 

0.22 

0 .IS 

0.22 

-0.12 

0.25 

0.26 

0.22 

0.26 . 

-0.06 

0.29 

0.30 

-0.16 

0.26 

-0.41 

0.32 

0.35 

-0.02 

0 . 32 

-0.39 

0.34 

0.57 

-0.27 

0.47 

-0.55 

O .44 

0.22 

-0.06 

0.23 

-0.20 

0 . 26 

0.19 

0.27 

0.20 

-0.16 

0.21 

0 .28 

0 . 08 

0.28 

-0.23 

0.31 

0.33 

-0.02 

0.34 

-0.06 

0 .32 

0.19 

-0.03 

0.18 

-0.43 

0.21 

0,18 

-0.33 

0.17 

-0.56 

0 .24 

0.20 

-0.19 

0.19 

-0.47 

0.25 


* Entire data set; rows not marked relate to the subset of data selected \o generate mode). 
} Solar Maximum Mission 
2 Channel ) 9 
: Channel 06 
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Table 7(b). Comparison of models to database on 0 7 concentration for high 
geomagnetic activity (A P > 50) 

Data Set ALT PTS MOO M90 J70 





ME AN 

St> 

MEAN 

3D 

MEAN 

SD 

SMK 1 

120 - 200 

4136* 

-0.25 

0.55 

-0,49 

0.61 

-0.69 

0.65 

UV Occ 

120 - 200 

1 S3 

-0.12 

0.54 

-0.44 

0 .46 

-0 . 67 

0.. 56 


200 - 400 

17 80* 

-0. B2 

3.77 

-0 . S3 

3 . S3 

-0.92 

3 .86 

AE-C OSS 

120 - 200 

2 

0.32 

0.23 

0.23 

0 .20 

0.42 

0.09 

ms 

200 - 400 

1 

1.28 . 

0 . 00 

1 . 00 

0.00 

0.70 

0.00 

AE-C EUVS 

100 

g 

0.06 

0.22 

0.06 

0.20 

-0.26 

0,15 

UV Occ 

130 

1 

0.26 

0.00 

0.19 

0.00 

0 . 00 

0 .00 


150 

7 

0.32 

0.39 

0.11 

0.39 

0,27 

0.38 

AS- 5 5UVS 

100 

3 

0 , 22 

0.19 

0.07 

0.14 

-0.05 

0.12 

UV Occ 

13 0. 

a 

0.09 

6.15 

0.04 

0 . 16 

-0.45 

0.21 


1 5 0 2 

5 

0 . 84 

1 . 61 

0.44 

1 . 67 

-0 . 19 

1 . 63 


150 3 

15 

-0.04 

0 .29 

-0.21 

0.30 

-0.53 

0.35 


? Emirs data set; rows noi marked relate lo the subs*! of data selected to generate model. 
J Solaj Maximum Mission 
7 Channel 3 9 
3 Channel 06 
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Table 7(c). Comparison of models no database on 0 2 concentration across all 
geomagnetic activity levels 

Data Set ALT PTS MOO M9 0 J7 0 





MEAN 

SD 

MEAN 

SD 

MEAN 

SD 

sm l 

120 - 200 

51 8 8 6 1 

-0.19 

0.37 

-0.42 

0.40 

-0.69 

0.44 

U V Occ 

120 - 200 

2091 

-0.09 

0,27 

» 

o 

UJ 

CD 

0.32 

-0.66 

0.39 


200 - 400 

27225- 

-0.25 

0.61 

-0.52 

0.72 

-0.72 

o.ao 

AE-C OSS 

120 - 2 0 D 

160 

0.07 

0.28 

-0.02 

0. 2B 

-0.13 

0.36 

N MS 

200 - 400 

66 

0.14 

0.39 

0 . 05 

0.37 

0.03 

0 .48 

AE-D OSS 

120 - 200 

376 

0.24 

0.34 

0.19 

0.34 

0.05 

0 .36 

N MS 

200 - 400 

133 

0.13 

0.32 

0.13 

0.33 

0.16 

0 . 45 

AE-S OSS 

120 - 200 

256 

0 . 2 B 

0.22 

0.22 

0.22 

-0.11 

0 . 25 

MM3 

200 - 400 

47 

0.15 

0.29 

0.15 

0.29 

-0.16 

0 . 32 

Rocket 

120 - 200 

14 

-0. D7 

0.30 

-0.12 

0.27 

-0.38 

0 . 32 

Absorption 









Rocket 

80 - 120 

37 

-0 . 17 

0.45 

-0,22 

0.45 

-0.44 

0 . 51 

MMS 

120 - 200 

295 

0.09 

0.34 

-0 . 04 

0.31 

-0.38 

0 . 33 


200 - 400 

44 

-0.19 

0.55 

-0.31 

0 . 45 

-0.56 

0 . 42 

AE-C EU VS 

100 

121 

0 . 00 

0.23 

0 . 04 

0.24 

-0.16 

0.25 

UV OCC 

130 

54 

0.28 

0.21 

0.26 

0.21 

-0.09 

0 . 24 


150 

130 

0 . 19 

0.31 

0.07 

0.31 

-0.15 

0.38 

AE-E EUVS 

1 0 D 

249 

0 . 06 

0.30 

0 . 03 

0 .31 

-0.05 

0 . 29 

UV Occ 

130 

175 

0 . 00 

0.20 

-0.03 

0.19 

-0.43 

0.22 


15 0* 

163 

-0.1S 

0.20 

-0.32 

0.19 

-0.56 

0 . 29 


150 3 

300 

-0 . 05 

0.21 

-0.13 

0.20 

-0.46 

0 .27 


; Entire daLa set; rows not marked relate to the subset of data selected to generate model. 

1 Solar Maximum Mission 

2 Channel 19 
? Channel 06 
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Table 8 (a) , 

, Comparison 

of models to 

database 

on atomic 

nitrogen 

concentration 


for quiet 

geomagnetic ac 

tivifcy levels (A p < 

; 10) * 


Data Sec 

ALT 

?TS 


1400 

M90 

J70 




MEAN 

SD 

MEAN 

SD 

MEAN SD 

AE-C OSS 

200 - 400 

380 

-0,11 

0.44 

0.11 

0.44 


NMS 

400 - 800 

1994 

-0 . 05 

0.30 

0 . 19 

0.31 


AE-D OSS 

200 - 400 

2160 

-0.42 

0.29 

-0.29 

0.30 


N MS 

400 - 800 

871 

-0.33 

0 ; 3 S 

-0.16 

0.36 


AE-E OSS 

200 - 400 

3998 

0.04 

0.27 

0.12 

0.27 


NMS 








DZ-2 NACS 

200 - 400 

169 

-0.75 

0.72 

-0.73 

0.68 


NHS 

400 - 800 

74 

-0.62 

0 .50 

-0.59 

0.47 



* Jacchia-70 model does not cover atomic nitrogen 
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Table 9(b). Comparison of models to database on atomic nitrogen concentration 
for high geomagnetic activity (h fi > 50}* 


Data 

Set 

ALT 

PT5 

MOO 

M30 

J7 0 





KEAN 

SD 

MEAN 

SD 

MEAN 

AE-C 

OSS 

200 - 400 

21 

-0 . 15 

0.36 

-0.02 

0.40 


N HS 


400 - 300 

335 

0. OC 

0 . 3 B 

0.19 

0.40 


DE-2 

NAC5 

200 - 400 

54 

-0.73 

1.13 

-0.69 

1.03 


NMS 


400 - 900 

26 

-0.07 

1.27 

-0.06 

1.19 



* Jacchia-70 model does not cover atomic nitrogen 
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Table S (c 


Data Set 


AE-C CSS 
NMS 


AE-D OSS 

ms 

AE-E OSS 
NMS 

DE-2 MACS 
NMS 


. ’ Jacchia 


Comparison 

of models to database 

on atomic 

nitrogen 

concentration 

across all 

geomagne 

tic scti 

.vity 

levels T 



ALT 

PTS 

,M0Q 

M90 

J7 0 



MEAN 

SD 

MEAN 

SD 

MEAN SD 

120 - 200 

2 

-0,96 

0.26 

-0.96 

0.29 


200 - 400 

670 

-0.12 

0.41 

0.09 

0.42 


400 - 800 

453 6 

-0 , 07. 

0.32 

0.17 

0.33 


200 - 400 

3 644 

-0.42 

0.31 

-0.28 

0.31 


400 - 800 

1659 

-0 . 37 

0.39 

-0.15 

0.39 


200 - 400 

7125 

0 . 16 

0 .30 

0.23 

0.30 


120 - 200 

103 

“0 . 44 

0.63 

-0.38 

0 . 63 


200 - 400 

903 

-0.61 

0.81 

-0 . 59 

0.73' 


400 - 800 

25 8 

-0:46 

0v7 3 

-O'. 4 4 

0.70 



70 model does not cover atomic nitrogen 
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Table 9(a) . 

Comparison 

of models to 

database 

on atomic 

hydrogen 

concentration 


for quiet 

geomagne 

tic activity 1 

evels (A ? < 

10) * 


Data Set 

ALT 

PTS 


MOO 

M90 

J70 




MEAN 

SD 

KEAN 

SD 

MEAN SD 

AZ-C BINS 

200 - 400 

1756 

0.03 

0.30 

0.05 

0 . 31 


IKS 

400 - 900 

466 

0.06 

0.40 

0.08 

0 . 41 


AZ-E BIMS 

120 - 200 

9 

2.7S 

0.11 

2.93 

0.12 


IKS 

200 - 400 

1113 

0.02 

0.29 

0.00 

0.27 



400 - 800 

516 

-0.09 

0.31 

-0 . 05 

0.32 



* Jacchia-70 model does not cover atomic hydrogen 


/ 
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Table 9(b). Comparison of models to database on atomic hydrogen concentration 
for high geomagnetic activity levels (A P £ 50) T 


Data 

Set 

ALT 

PTS 

H00 


M9 0 





MEAN 

SD 

MEAN 

SD 

AE-C 

3IKS 

200 - 400 

220 

0.04 

0.59 

0.04 

0. 60 

IMS 


400 - 800 

S 5 

-0 . 01 

0.37 

-0.02 

0.36 

AE-S 

BIMS 

200 - 400 

71 

0 .16 

0.40 

0.13 

0.47 

IMS 


400 - 800 

33 

-0.19 

0.29 

-0.15 

0.29 


* Jacchia-70 model does not cover atomic hydrogen 


J7 0 

MIAN SD 
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Table 9(c). Comparison of models to database on anomic hydrogen concentration 
across all geomagnetic activity levels* 


Data 

Set 

ALT 

?TS 

A3-C 

3IMS 

200 - 400 

3735 

IMS 


400 - BOO 

942 

AS-5 

3IMS 

120 - 200 

16 

IMS 


200 - 400 

22^5 



400 - 800 

1084 


MOD M90 


MEAN 

SD 

MEAN 

SD 

0.06 

0.32 

0.07 

0.33 

0 . 03 

0 . 3 B 

0 . D3 

0 . 3 B 

2 .10 

0.31 

2.24 

0.31 

0 . 01 

0.30 

-0.01 

0.30 

-0. 11 

0.34 

-0.07 

0.35 


J7 0 

MEAN SD 


* Jacchia-70 model does not cover atomic hydrogen 
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POPULAR SUMMARY 


A new global, upper-atmospheric model is presented. This model is intended to replace 
previous models utilized in such applications as satellite orbit prediction. The new model 
termed the NRLMSISE-00 model, together with an associated database, NRLMSIS, is 
based on total mass densities from satellite accelerometers, satellite orbits, temperature 
determined by incoherent scatter, and molecular oxygen densities from the SMM solar 
occultation experiment. Extensive comparison is made with earlier upper atmospheric 
models such as MSIS and Jacchia-70. Model dependence on solar and geomagnetic 
activity is demonstrated more clearly than previous models. This model should become 
the standard for use in satellite orbit prediction. 



